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TWO-DIMENSIONAL  TRANSONIC  FLOW  PATTERNS  *  1 

By  Stefan  Bergman. 

1.  Introduction.  The  successful  development  of  the  mathematical 
theory  of  two-dimensional  steady  flows  of  an  ideal  incompressible  fluid  is 
largely  due  to  the  fact  that  the  complex  potentials  of  such  flows  are  analytic 
functions  of  a  complex  variable.  Function-theoretical  methods  may  therefore 
be  used  in  order  to  obtain  and  to  investigate  flow  patterns  possessing  at  given 
points  the  assigned  singular  character  and  satisfying  given  boundary 
conditions.2 

In  investigating  the  flow  of  an  incompressible  fluid  by  theoretical  methods, 
two  alternative  treatments  have  been  used.  The  behavior  of  the  solutions  has 
been  studied  either  in  the  so-called  physical,  i.  e.,  the  plane  in  which  the  motion 
actually  occurs,  or  in  the  hodograph  plane,  i.  e.,  the  plane  whose  cartesian 
coordinates  are  the  velocity  components.  In  the  incompressible  case,  both  the 
'potential  and  the  streamf unction  are  harmonic  functions,  irrespective  of 
whether  the  motion  is  considered  in  the  physical — or  the  hodograph  plane. 

In  .  the  compressible  case,  both  the  potential  and  the  streamfunction 
considered  in  the  physical  plane  satisfy  complicated  non-linear  equations.  By 
considering  the  motion  in  the  hodograph  plane,  and  making  a  few  appropriate 
transformations,  it  is  possible  to  linearize  these  equations.  (See  [11]  and  2 
of  the  present  paper.3) 

The  treatment  of  the  problem  in  the  physical  plane  has,  of  course,  the 
fundamental  advantage  of  making  the  boundary  conditions  enter  in  an  obvious 
way.  The  gain  in  simplicity  due  to  the  linearization  is,  however,  so  consider¬ 
able  as  to  make  the  hodograph  method  vastly  superior  for  various  purposes. 

In  the  present  paper  the  theory  of  compressible  fluid  flow  is  developed 

*  Received  August  2,  1947. 

1  Research  paper  done  under  Navy  Contract  NOrd  8555-Task  F,  at  Harvard  Uni¬ 
versity.  The  ideas  expressed  in  this  paper  represent  the  personal  views  of  the  author, 
and  are  not  necessarily  those  of  the  Bureau  of  Ordnance. 

2  Poles,  logarithmic  singularities,  etc.,  of  the  complex  potential  represent  doublets, 
sinks,  sources,  vortices,  etc.,  of  the  flow. 

3  The  numbers  in  brackets  refer  to  the  bibliography  at  the  end  of  the  paper. 
Acquaintance  with  the  contents  of  these  publications  is  not  assumed  in  the  present 
paper. 
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by  studying  the  linearized  equation  and  applying  to  it  the  operator  method 
developed  in  [1-10].  Generalizing  the  procedure:  “taking  the  real  part/5 
this  method  leads  to  the  determination  of  certain  linear  operators  which  trans¬ 
form  functions  of  one  variable  into  stream  (or  potential)  functions  of  com¬ 
pressible  fluid  flows  preserving  many  fundamental  properties  of  functions  to 
which  the  operator  is  applied.  In  the  present  paper,  two  operators  of  this 
kind  are  considered:  the  so-called  integral  operator  of  the  first  kind  (3)  and 
that  of  the  second  kind  (4-6).  For  some  fluid  dynamical  applications,  and 
in  particular  for  the  study  of  the  behavior  of  the  flow  near  the  sonic  line, 
the  second  operator  seems  to  be  more  appropriate.  In  this  case,  however,  the 
relations  between  the  properties  of  functions  to  which  the  operator  is  applied 
and  those  of  generated  functions  is  more  hidden.  This  is  why  it  is  useful  at 
first  to  investigate  the  special  case  which  is  obtained  by  assuming  a  simplified 
equation  of  state.  This  simplifying  assumption  results  in  the  equation  for 
the  stream  function  taking  the  form 

(1*1)  —  0H(0V/M)  +  (0V/0H2)  -0,  C>  0 

where  H  is  a  function  of  the  Mach  number,  M ,  which  is  negative  for  M  <  1 
and  positive  for  M  >  1.  (See  (2. 1)  and  (2.  3)). 

In  an  important  investigation,  Tricomi  [16]  studied  the  boundary  value 
problem  of  equation  (1.1)  and  showed  that  if  we  consider  a  finite  domain  D, 
bounded  in  the  supersonic  region  by  two  characteristics,  say  BA  and  GA , 
and  by  a  curve  BmO  in  the  subsonic  region,  and  if  the  boundary  values  are 
prescribed  in  BmC  and  on  one  of  the  characteristics,  say  BA,  then  the  boundary 
value  problem  has  a  unique  solution.  Frankl  [13]  considered  questions  allied 
with  Tricomfis  investigations  in  the  case  of  the  exact  compressibility  equation. 

The  questions  arising  in  our  approach  are,  however,  of  a  somewhat 
different  nature  than  those  considered  by  the  above-mentioned  authors. 

In  the  first  place,  we  seek  to  find  conditions  for  a  function,  say  /,  of  one 
variable  in  order  that  the  generated  function  P(/)  will  be  defined  in  a  pre¬ 
scribed  domain  B ,  which  in  general  lies  partially  in  the  subsonic  and  partially 
in  the  supersonic  region. 

Secondly,  and  this  is  the  most  essential  difference,  we  are  considering 
solutions  of  the  compressibility  equation  which  possess  singularities  (e.  g., 
branchpoints)  in  the  hodograph  plane.  In  the  applications  of  the  theory, 
the  consideration  of  this  type  of  solution  cannot  be  dispensed  with,  since, 
notwithstanding  the  singularities  in  the  hodograph  plane,  the  behavior  of  the 
solutions  in  the  physical  plane  can  be  perfectly  regular.  Furthermore,  certain 
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singularities  in  the  physical  plane  have  a  hydrodynamical  meaning  and  must 
be  considered  in  investigating  flow  patterns. 

The  methods  employed  in  the  case  of  equation  (1.1)  are  to  some  extent 
capable  of  generalization  to  the  case  in  which  the  coefficient  ( —  (7H)  is 
replaced  by  an  arbitrary  function  l (H) .  (See  5) .  This  includes,  in  particular, 
the  exact,  i.  e.,  non-simplified  compressibility  equation. 

In  3,  we  introduce  the  so-called  integral  operator  of  the  first  kind,  Pi. 
This  operator  yields  a  streamfunction  of  a  subsonic,  compressible  fluid  flow 
in  terms  of  an  arbitrary  function  of  one  complex  variable.  The  representation 
holds  for  4  E[M  <  1, —  <  0  <  oo].  Here  M  is  the  Mach  number  and  0 

the  angle  which  the  velocity  vector  forms  with  the  positive  a-axis. 

An  analogous  representation  for  the  streamfunctions  of  supersonic  flows 
in  terms  of  two  differentiable  functions  of  one  real  variable  holds  for 
E  [M  >  1, —  oo  <  0  <  oo  ].  Using  the  integral  operator  of  the  second  kind, 
we  obtain  (4,  5)  four  analogous  representations  in  terms  of  arbitrary  functions 
of  one  variable.  These  four  representations  are  valid  in  four  adjacent  domains 
of  the  M,  0-plane,  namely 

D1  =  E[M  <  1,0  >  34  |  A (M)  |]  +  E [M  >  1,  6  >  A (M)], 

B2  =  E [M  <  1,  |  0  |  <  3^  |  \(M)  |], 

Z>3  =  E[M<  1,  0  <  —  |  X(M)  |]  +E \_M  >  1,  0  <  3A(Af)  ], 

D4  =  E[Af  >  1,  — 3A (M)  <  0<  A(ilf)], 

respectively.  Here  0=±3^A(Af)  and  0=  ( — lzk2)A(M)  are  certain 
curves  which  pass  through  the  point  M  =  1,  0=0  and  which  lie  in  the  sub¬ 
sonic  and  supersonic  region  respectively  (see  fig.  1). 

In  the  simplified  case  these  four  representations  can  be  combined  into 
one,  yielding  a  representation  which  holds  in  the  whole  M ,  0-plane.  This 
result  is  based  on  certain  theorems  of  the  Fuchs  theory  of  ordinary  differential 
equations  with  singular  coefficients.  The  question  of  combining  the  analogues 
of  the  above  four  solutions  in  the  exact  case,  and  generally,  the  study  of  the 
solutions  leads  to  the  investigation  of  partial  differential  equations  with 
singular  coefficients,  which,  when  solutions  are  continued  to  complex  values 
of  the  arguments,  can  be  attacked  by  methods  representing  a  generalization 

4  The  functions  t  =  PA/)  may  he  multi-valued  functions  which  may  possess  singu¬ 
larities.  The  statement  that  P Af)  is  defined  in  E[ilf  <  1]  means  that  the  projections 
of  the  domain  in  which  P 1{f)  is  defined  on  the  schlicht  M,  (9-plane,  lies  in  E[3f  <1]. 

E[  ]  denotes  the  set  of  points  whose  coordinates  satisfy  .conditions  indicated  in 
the  brackets. 
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of  the  Fuchs  theory  for  ordinary  differential  equations.  These  questions,  and 
in  particular  the  problem  of  combining  the  four  above  representations  into 
one,  will  be  treated  in  a  subsequent  paper. 

In  6,  we  determine  the  “associate”  function  for  P2(/)  in  terms  of  the 
values  of  the  stream  function  ^  =  Im[P2(/)],  and  its  derivative  with  respect 
to  M  on  the  line  M  =  1  (sonic  line). 

The  author  would  like  to  take  this  opportunity  to  thank  Bernard  Epstein 
and  A.  Zeichner  for  helpful  advice  and  aid  in  connection  with  the  present 


Fig.  I. 


paper.  He  would  also  like  to  thank  Z.  ISTehari  and  M.  Schiffer  for  a  number 
of  helpful  suggestions. 

2.  Equation  for  the  streamfunction  of  a  compressible  fluid  flow. 
Exact  and  simplified  equations.  Assuming  that  the  thermodynamical  equa¬ 
tion  of  state  of  the  fluid  has  the  form  p  —  o-pk ,  where  a  and  Tc  are  constants, 
and  p  and  p  the  density  and  pressure  respectively,  and  introducing  as  new 
variables 

(2.1)  H=  rP(dq/q)=  ff[l  -  J(*  —  l)?2]1^1^ 

and  6,  where  q  is  the  speed  and  0  the  angle  which  the  velocity  vector  makes 
with  some  fixed  direction  (say  the  positive  ®-axis,),  we  obtain  the  following 
linear  equation  for  the  streamfunction : 
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'  (2,2)  S(t)  =Z(H)(0V/^2)  +  mm2)  -0,  1(H)  =  (1—  M2)/p2, 

where 

(2.  3)  M  =  q/[l  —  i(k  —  l)q2]h 

is  the  Mach  number.  The  denominator  in  (2.  3),  being  the  local  velocity  of 
sound,  M  will  be  smaller  or  larger  than  1  according  as  the  flow  is  subsonic 
or  supersonic  respectively.  The  differential  equation  (2.  2)  will  therefore 
be  of  ellipitic  or  hyperbolic  type,  corresponding  to  the  subsonic  or  supersonic 
character  of  the  flow. 

A  formal  computation5  shows  that  the  Taylor  development  of  l( H)  in 
the  neighborhood  of  H  —  0  is 

(2.4)  Z(H)  —  [0/(fc  —  —  2H) 

—  ( (21c  +  5) /(2k  +  2))((k  +  1)  /2)2k/ik~1)  (—  2H)2 

_j_  A;4  -[-  (43/6) k3  -f~  -j~  (31/2) fc  fl-  31/6  _|_  q) /2)fc/(ft_1)  ( 2H)3  -|- 

(k  +  l)4 

In  considering  a  flow  (or  a  portion  of  a  flow)  which  is  purely  subsonic, 
it  has  some  advantages  to  replace  H  by  the  variable  A  defined  by 

r_H 

(2.5)  — A—  'I  P(—T)]»tfr. 

r= o 


A  can  be  expressed  in  a  closed  form  as  a  function  of  M;  a  formal  compu¬ 
tation  yields 

.  r*i —  (i-^-M2)*  (i  +  h(i—M2)*Yk-\ 

(3.6)  .  A-  2  log  Li  +  (1  —  (l  —  h(l—  M*)i)  J’ 

h  =  [(h- l)/(i  +  l)]». 

Now 

—  i/oAh  *=  —  ^xpZ(H)]^, 

(2. 7)  ^Hh=  «AxxZ  +  ^x(^)  h  =  '/'xx*  —  i^xHZH  =  Zfyxx  +  ^"3/2^h  ] 


so  that  (2.  2)  becomes 

(2.  8a)  ^xx  +  \jjQQ  +  =  0, 


e  A  detailed  account  of  formal  derivations  of  some  expressions  used  in  the  present 
paper  can  he  found  in  the  Appendix  to  Technical  Report  10,  of  the  series  “  Operator 
methods  in  the  theory  of  compressible  fluids,”  Harvard  University,  1948. 
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or,  in  complex  notation, 

(2-  8b)  fyzz  +  4 N (if/z  +  $z)  =  0, 

where 

(3.9)  ff_jMH__[(i  +  1)/g]__^_i  + 

See  [6,  (46)].  It  should  be  noted  that  the  interval  —  oo  <  \  <  0  corresponds 
to  the  interval  —  oo  <  H  <  0. 

In  the  supersonic  case  (i.  e.,  for  M  >  1)  the  right-hand  side  of  (2.  6) 
becomes  purely  imaginary.  If  we  introduce  the  new  variable  A  defined  by 

(2.10)  A  — 
it  is  easily  confirmed  that 

(2. 11)  A  =  hr1  arctan  [h(M2  —  1)3]  —  arctan  [( M 2  —  1)*]. 

In  this  case,  (2.  8a)  will  take  the  form 

(3.13)  *AA-^  +  4^A=0, 

Remark.  Equations  (3.8)  and  (3.13)  can  be  simplified.  If  ^  is 
replaced  by 

(3. 13)  =  MR 

where  ( SR/dZ )  —N,  then  ifi*  satisfies  the  equations 

(3. 14)  +  ^ee  +  4Fty*  =  0  and  —  4 Ftf*  -  0, 

where 

(2  ik\  R  =  R  __  f)4f4  p  (3fc  1  )4f4  —  4(3  —  3fc)itf2  -j-  16"] 

'  1  64  L  (1  —M2Y  J‘ 

For  subsequent  use,  we  write  down  the  expansions  of  N  and  F  in  the 
neighborhood  of  A  =  0 

(3.16)  tf=(l/m)[l  —  |(J:  +  l)![(3  +  -|3S)fc+5-3i  —  3]  X 

[3-l/«32/3(J.  _|_  !)-./.(_  X)2/3] 

(3.17)  F=  (5/144)  ( —  A,)'2  +  4.2( — A)_2/s  +  J.0  +  A2(- — ■A)2''3 

+  •  •  •  =  5/36  ( —  3A)2  +  ' 

In  the  vicinity  of  H  =  0,  i.  e.,  the  sonic  line,  1(H)  may  be  replaced  by 
the  first  term  in  its  expansion  (3.4).  Using  this  value  of  1(H)  in  (2.3), 
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we  obtain  the  so-called  “simplified”  compressibility  equation  (1.1).  In 
considering  transonic  flows,  the  solutions  of  the  simplified  equation  will 
therefore  give  a  fair  approximation — in  a  certain  neighborhood  of  the  sonic 
line— of  the  exact  streamfunction. 

The  expression  for  N,  F,  and  H  will,  in  this  case,  reduce  to 

(2.18)  N  —  Wf  —  —  (1/6)  (1/ ( —  2A) ), 

P  =  Pf  -  (5/36)  (l/( —  2A)2), 

H  =  Hf  =  (32/3/2)  (2/ (7c  —  1)  )  (fc-2>/(3fc-3)  (—A)2/3, 

respectively. 

Kemark.  We  are  using  the  same  variable  A  in  both  the  exact  and  the 
simplified  case,  as  this  facilitates  the  comparison  of  the  respective  flow  patterns. 

3.  Application  of  integral  operators  to  the  compressibility  equation. 
Integral  operator  of  the  first  kind.  The  use  of  integral  operators  in  the 
theory  of  the  compressibility  equation  is  based  on  the  following  theorem: 

Theorem  3. 1.  Let  E(Z,Z,t)  be  a  function  of  two  real  and  one  com¬ 
plex  variables ,  A,  6,  t ,  which  is  defined  for  t  along  a  curve  connecting  t  —■  1 

and  t  =  i9  and  for  (a,  6)  e  G.  G  denotes  here  a  sufficiently  small  neighborhood 
of  the  origin. 

Let  E  satisfy  the  following  conditions: 

1.  E  possesses  continuous  partial  derivatives  with  respect  to  all  three 
of  its  arguments ,  up  to  the  second  order. 

2.  The  expression 

(3.1)  [(i  —  P)E{Z,z9t)/Zf\d/M 

is  continuous  for  Z  =  0,  and  approaches  zero ,  uniformly  with  respect  to 
(A,  6)  £  G  as  t  -»  —  1  or  f  -»  1. 

3.  E  satisfies  the  partial  differential  equation 

(3.  2)  G(E)  as  (1  —  t 2)  (Bit  +  NEt)  -  (1  /t)Ei  +  2 tZL(E)  =  0 
where 

(3.  3)  L(E)  aa EZz  +  N(EZ  +  El)  as  (l/4)2?xx  +  (1/4)  +  NEX. 

If  f(£/2)  is  an  analytic  function  of  £  defined  in  a  simply-connected  domain  P, 
which  includes  the  origin ,  then  the  expression  u(Z,Z),  given  by 
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(3.4a)  u{Z,Z)=V{f), 

(3.  4b)  P  (/)  =  j\E(z>  Z>  *)f(iz(  1  -  n  )<«/(!  -  *2)s 

defined  in  a  simply-connected  domain  which  lies  in  G  f)  P,  and  satisfies 
the  equation  L(u)  =0. 

The  proof  of  this  theorem  is  given  in  [1,  §  1  and  6,  pp.  34-39].  The 
function  /  will  be  called  the  associate  of  P(/)  with  regard  to  the  operator  P. 

Since  (3.  2)  has  an  infinity  of  solutions,  there  exist  infinitely  many 
integral  operators;  a  closer  investigation  of  their  properties  will  show  which 
type  of 'operator  is  best  suited  for  the  purpose  on  hand.  The  following 
property  of  the  integral  operator  plays  an  important  role  in  some  of  the 
applications. 

As  is  well-known,  a  harmonic  function  Gx( X'fl)  =  G(Z',Z')  may  be 
written  in  the  form 

G(Z',Z')  =  (l/2i)[g(Z')-g(Z')]. 

Here  6 

Z'  —  iO' ,  ZJ  =  X'  —  i0' ,  X.'  —  A  —  A0,  0'  —  0  — -  60 ,  A0  <C  0. 

If  now  the  harmonic  function  Gt( A',  O')  is  continued  to  the  complex  values 
of  the  arguments  A'  and  O',  i.  e.,  if  we  assume  that  Z'  and  Z'  are  not  necessarily 
conjugate  to  each  other,  and  if,  in  particular,  we  consider  G  and  g  in  the 
so-called  characteristic  planes  Zr  =  0  or  Zf  =  0,  then  we  see  that  the  analytic 
function  of  a  complex  variable  and  the  continuation  of  the  real  harmonic 
function  differ  only  by  constants ; 7  indeed, 

G(Z',0)  =  (l/2i)[g(Z')-g(0)l 

An  integral  operator  generates  complex  solutions  of  L,  and  we  may  demand 
that  these  complex  solutions  possess  an  analogous  property.  We  shall  show 
that  there  exists  an  integral  operator  p  (which  is  connected  by  relations 
(3.  11),  (3. 12)  with  P)  such  that  the  complex  solution  of  (3.  3) 

0  In  the  discussions,  it  will  he  useful  to  consider  a  shift  of  origin  to  the  point  \0,  0O 
(point  of  reference  of  the  operator). 

In  this  section  Z'  and  Z'  will  be  treated  as  independent  variables. 

7  Analytic  functions  of  a  complex  variable  represent  a  very  special  subclass  of 
complex  harmonic  functions  (i.  e.,  the  totality  of  functions  G  +  ill ,  where  G  and  H  are 
two  arbitrary  real  harmonic  functions).  Operator  (3.4a),  (3.4b)  generates  also  a 
subclass  of  complex  solutions  of  the  equation  £. 
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u{Z',Z')=pig{Z')\ 

and  the  real  solution 

t(Z',Z')  =  lm[u{Z',Z/)'] 
are  connected  by  the  following  relations 

(3.  5)  $(Z',  0)  =  (1  J2i)  [u (Zf,  0)  —  R{Z',  0)  const.],  u(Z'} 0)  =  g(Z') 

where  R(Z',Z' )  is  a  given  function  defined  in  (3.  7).  This  operator  P  will 
be  called  “integral  operator  of  the  first  kind”  and  will  be  denoted  by  Pi. 

Defining  the  generating  function  Ex{Z\rZ/,t)  of  the  operator  of  the 
first  kind  by  the  requirement  that 

(3.6a)  Ei(Z',0,t)  =1 

and 

(3.  6b)  ^(O,  Z',  t)  =  exp  [  —  C  NdZ% 

J  0 

we  shall  show  that  these  relations  imply  the  property  (3.  5). 

Writing  the  generating  function  of  the  first  kind  E±  in  the  form 
(3.7)  E^RiZ'J^E^Z'J'J); 

R(Z',Z')  =  exp  [—  f*N(Z'  +  Z'1)dZ'1] 
and  assuming  that  E*\  has  the  development 
(3.  8)  E*1  =  1  +  1  {Z\  Z'), 

n= l 

it  is  found  that  (3.7)  satisfies  the  relation  (3.6b).  Substituting  this  into 
equation  (3.  2)  we  find  that  the  P(n)  satisfy  the  following  recurrence  relations : 

(3.  9)  Pz>(1)  +  2F  =  0,  (2n  -j-  1  )P2'(n+1)  +  2PZ’Z'{n)  +  2 F =  0, 

ft  =  1,  2,  3,  •  • 

Finally,  (3.  6a)  is  satisfied  by  imposing  upon  the  P{n)  the  initial  conditions 

(3.  10)  P<»>  ( Z',  0)  =  0,  ft  =  1,  2,  3,  ■  ■  • . 

By  the  above  requirement,  the  P{n)  and  hence  the  generating  functions  E ± 
(of  the  first  kind),  are  uniquely  determined.  Applying  the  considerations 
of  [1,  pp.  1173-76]  it  can  be  shown  that  the  series  (3.  8)  converges  absolutely 
and  uniformly  in  a  sufficiently  small  neighborhood  of  the  origin  Z/  =  0, 


TWO-DIMENSIONAL  TRANSONIC  ELOW  PATTERNS. 


865 


Z'  =  0.  The  existence  of  an  integral  operator  of  the  first  kind  thus  is 
assured.  Assuming  that  the  associate  function  /  is  regular  in  a  sufficiently 
large  domain,  we  prove  that  by  applying  to  it  the  integral  operator  of  the 
first  kind,  we  obtain  a  solution  u(Z,Z)  (see  (3.4a),  (3.4b))  of  (2.8b) 
defined  in  a  sufficiently  small  neighborhood  of  the  origin.  We  shall  show  in 
the  following  that  if  /  is  regular  for  M  <  1,  this  solution  can  be  continued 
throughout  the  whole  subsonic  region. 

Remark.  Integral  operators  of  the  first  kind  can  also  be  written  in  a 
somewhat  different  form  which  is  useful  for  various  purposes.  Namely  (as 
can  be  shown  by  a  straightforward  computation,  see  [4,  pp.  618-619]) 
we  have  8 

(3.11)  p  (f)-p(g) 

=  R(Z',Z')[g(Z')  +  22-*»  rr(,2”+ V  PM(Z',Z')gM(Z')] 

n= l  1  \n  l  *-) 

Z'  />  Z i  Zn-i 

(  g{Zn)dZn •  •  ■dZ1 

0  t/  0  0 

(-l)n- 1  r  r=Z' 

“fc=ip  J* 

where 

(3.12)  <?(Z')  =  r,/[P,(l-<2)]^/(l-i2)s- 

*/  f=-l 

(3.11),  (3.7)  and  (3.10)  imply  (3.5). 

We  proceed  now  to  the  proof  that  every  real  solution  of  equation  (2.  8b) 
can  be  represented  in  a  sufficiently  small  neighborhood  of  the  origin  Z '  =  0, 
Z'=0  as  the  imaginary  part  of  the  right-hand  side  of  (3. 11)  with  suitably 
chosen  associate  function  /  (or  g). 

Let  \p(Z',Z')  (ZZ  being  conjugate  to  Z')  be  a  real  solution  of  equation 
(2.  8b)  which  is  regular  in  a  sufficiently  small  neighborhood  of  the  origin. 
Since  this  equation  is  of  elliptic  type  and  its  coefficient  N  is  an  analytic 
function  of  two  variables,  if/(Z',Z')  can  be  written  in  the*form  of  a  power 
series 

00  00 

<P{Z',  Z')  =  2  2  DmnZ'mZ'n,  Dmn  =  Dnm, 

m= o  n- o 

which  converges  in  a  sufficiently  small  neighborhood  of  the  origin. 

8  f  and  g  are  associates  of  the  same  solution  of  (  2.  8b),  the  first  with  respect  to  the 
operator  P,  the  second  with  respect  to  p.  In  order  to  avoid  any  confusion  we  shall 
speak  about  “  p-associate  ”  and  “  P-associate,” 
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G1  and  G2  are  two  analytic  functions  of  one  complex  variable  Z'  and  Z\ 
respectively,  which  are  regular  in  a  sufficiently  small  neighborhood  of  the 
origin.  (We  note  that  Gi(0)  =Cr2(0)  and  G2(Z')  —  Gi(Zf)  since  for  rZ' 
conjugate  to  Zr,  ij/  is  real.)  On  the  other  hand,  by  classical  results  (the  initial 
value  problem  in  the  theory  of  partial  differential  equations),  it  is  known 
that  if  functions  G1(Z/);>  G2(Z'),  Cri(O)  =  G2(0)y  are  given,  there  exists  one 
and  only  one  solution  \\i{Z',Zr)  of  equation  (2.8b)  such  that  (3.13)  and 

(3.14)  hold.  The  integral  operator  (3.4)  enables  us  to  write  down  the 
solution.  Indeed,  let  us  determine  two  functions  say  gi(Z')  and  gz{Z'), 
g2(Zf)  =^g1(Z/),  such  that 

(3. 15)  9l(Z')  +  g2( 0)WJ')  =  W) 
and  therefore : 

g2(Z')  +  #i(0)I?(0,  Zf)  =  G2(Z') . 

Now 


(3.16)  R(Z',Z')\gi(Z') 


<~<*»>  r  ■  ■  r 


+  R{Z',Z')\g2{Z') 


,  £  T(2n  +  1) 
^*2*»T(n  +  1) 


will  represent  a  solution  of  (2.  8b)  which  satisfies  the  conditions  (3. 13)  and 
(3. 14)  ;  our  assertion  that  every  (real)  solution  can  be  represented  as  the 
imaginary  part  of  (3. 13)  is  therefore  proved. 

As  already  mentioned,  (3. 18)  is  a  priori  only  defined  in  a  sufficiently 
small  neighborhood  of  the  origin.  We  shall  prove,  however,  that  provided  f 
is  regular  for  A' <  A0j  the  solution  \f,(Z\Z')  obtained  in  this  way  can  be 
continued  into  the  whole  region  E  [Re(Z7  -J-  Zf)  <  2A0,  |  Zf  \  <  co,  |  Zf  |  <  oo]. 
As  shown  in  [6,  pp.  56  ff. ;  8,  p.  48  ff.],  the  quantity  F  =  F(^(Z'  -f-  Z') ) 


TWO-DIMENSIONAL  TRANSONIC  PLOW  PATTERNS. 


867 


introduced  in  (2. 15)  is  a  function  of  two  complex  variables  Z\  Zr  which  is 
defined  in  the  above  region.  Therefore  the  expression  9 

(3.17)  R(Z',Z')[gi(Z')-~  fZ  fZFg1dZ1dZ1 

J  O  J  0 

+  fZ  fZl  (  ZFg1dZ2dZ2]dZ1dZ1  +  -  ■  .] 

«/  0  ^  0  t/  o  ^  0 

+  B(Z',  Z')  [g2(Z ')  -  C’  f  Fg2dZ\dZ\ 

J  0  J  o 

H~  f*  (*  (*  r~  Fg2dZ2dZ2]dZ1dZ1  -j-  •  •  •] 

Jo  J  o  J  o  J  o 

satisfies  the  differential  equation  (2.  8b)  and  the  initial  conditions  (3. 13), 
(3.14).  It  is  evident  that  the  series  (3.17)  converges  in  any  'simply- 
connected  domain'  which  includes  the  origin  Z'  =  0,  Z'  =  0  and  which  is 
common  to  the  regularity  domains  of  F,  gx,  g2.  Since  by  the  above  require¬ 
ments  a  solution  of  (2.  8b)  is  uniquely  determined,  the  expressions  (3. 16) 
and  (3. 17)  must  coincide,  so  that  they  are  two  different  representations  of 
the  same  function. 

We  proceed  now  to  the  discussion  of  the  relations  between  the  domains 
of  regularity,  a  and  Tc ,  of  the  pr -associate  function  g  and  the  generated 
solution  pi(g)  in  the  real  plane,  i.  e.,  for  Zr  conjugate  to  7/. 

Theorem  3.  2.  Let  B  be  a  hounded  region  of  the  (real)  A',  6' -plane, 
situated  in  E[A'  <  A0].  If  g  is  regular  in  B ,  then  px(g)  is  also  regular  in  B ; 
conversely,  the  regularity  of  px(g)  in  B  implies  that  g  is  regular  there. 

Proof.  In  order  to  prove  our  statement  we  investigate  the  relations 
which  exist  between  the  regularity  domain  of  a  solution  \\i(Zf,Z')  of  (2.  8b) 
in  the  real  A',  0'-plane  (i.  e.,  for  Zf  conjugate  to  Z')  and  in  the  space  of  two 
complex  variables,  A,  6  ( i.  e.,  when  Z'  and  Zf  are  two  independent  complex 
variables).  Let  B  be  a  domain  in  the  A,  0-plane.  We  denote  as  the  hull10 

9  The  integral  operator  of  the  first  kind  may  be  regarded  as  a  generalization  of  the 
Riemann  formula  in  the  theory  of  linear  hyperbolic  equations  uXy  +  »(X,  Y)ux 
+  b  {X,  Y)  uY  +  c  (X,  Y)  —  0  to  the  elliptic  ease,  where  the  real  variables  X,  Y  are 
replaced  by  two  independent  complex  variables,  Z'  and  Z'  respectively.  If  we  consider 
the  solution  in  the  real  plane,  i.  e.,  for  Zf,  //  conjugate,  we  thus  obtain  solutions  of 
elliptic  equations.  See  [5,  pp.  317-318]. 

10  The  superscript  indicates  the  dimension  of  the  manifold  under  consideration. 
In  the  case  where  the  manifolds  are  one-  or  two-dimensional  and  are  situated  in  the 
(real)  X,  0-plane,  these  superscripts  are  omitted. 
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H4(B)  of  B  the  intersection  B x4  f]  B2 4  of  two  cylinders  of  (fonr-dimensional) 
space,  B14=E[Z'  e  B,  2'  arbitrary],  £24  =  B[Z/  arbitrary,  Z'eB]. 

It  is  well-known  that  a  solution  Zf)  of  the  first  equation  of  (2. 14) 

can  be  represented  in  the  domain  B  of  the  real  A,  0-plane  in  the  form 

(3. 18)  r{Z', z')  =  £ z'l LX)  —r(L Z)  g)-] 

where  b  denotes  the  boundary  of  B ,  the  interior  normal  to  b,  ds$  the  line 
element  of  b.  Here  <£*  is  a  fundamental  solution  of  (2. 14),  and  we  have 

(3.19)  <*.*  =  4x(Z',z';£,?)  [lg(Z'—0  +  lg(z'  — 
where  11 

(3.20) -  X(Z',Z';U) 


rz'  r2' 

—  1  —  j  j  FdZxdZj 

'  +  MW 


FdZ2dz£\  dZxdZx  -}-  * 


(3.21)  v(Z',Z'U,l) 


J'Z'  rz* 

I  GdZ1dz1 

r 

-srsr^cs: 


GdZ2dz£\  dZxdZx  -J- 


G  =  -(l/(Z'-l))  (dx/dZ')  -  ( 1  /(Z'-Q)  (dx/dZ') 

Since  F  is  defined  for  all  values :  Ee((Z/  +  Z')/2)  <  A0,  |  Z  |  <  co,  |  Z  |  <  co; 
it  follows  from  (3. 17)  that  if  gx  is  regular  in  B,  B  =  E[A'  <  A0],  then  pi(gi) 
is  regular  in  E4(B),  and  therefore  also  in  the  domain  B,  which  is  the  inter¬ 
section  of  H4(B)  with  the  real  A,  0-plane.  This  is  the  first  assertion  of 
Theorem  3.  2. 

From  (3.18),  (3.19),  (3.20),  (3.21)  it  follows  that  every  solution 
which  is  regular  in  the  domain  B  of  the  real  A,  0-plane  can  be  extended  to 
the  complex  values  in  H4(B).  Since  B  represents  the  intersection  of  H4(B) 
with  the  plane  =  0  as  well  as  with  Z' =  0,  ij/(Z',  0)  and  i[/(0,Z')  are 
regular  in  B.  Since  R(Z/,  Z')  is  regular  in  B[Re((Z/  -j-  Zf)/2)  <  A0,  |  Zr  |  <  oo, 
|  Zr  |  <  oo ],  it  follows  from  (3. 17)  that  Gx{Zf),  (as  well  as  Gz{Z')  =  GX(Z') 
see  (3.15))  is  regular  in  B.  See  [1,  §1  and  2,  §2].  This  completes  the 
proof  of  Theorem  3.  2. 

11  We  note  that  0,0)  =  />,(!). 
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As  has  been  shown  in  [5,  pp.  318-319],  the  function  generated  by  integral 
operators  P  (/)  defined  in  (3. 11),  will  have  branchpoints  of  finite  order,12 
at  the  same  points  as  the  function 

(3.22)  Po(/)=  fl/[iZ(l  —  t*)]dt/(l  —  t*)K 

Jt=- 1 

As  a  consequence,  the  following  holds. 

Theorem  3.  3.  Suppose  that  the  function 

(3.23)  9(Z)-VQ(f) 

is  defined  and  regular  in  a  region  R  situated  on  a  Riemann  surface ,  which 
possesses  in  its  interior  a  finite  number  of  branchpoints ,  each  of  finite  order . 
Let  further  the  projection  of  R  on  the  schlicht  A,  6-plane  lie  in  E[ —  oo  <  A' 
<A„]. 

Then  the  function 

(3.24)  Pi(/)=  — 

J  t=- 1 

(see  (3.2),  (3.7),  (3.8),  (3.9),  (3.10)),  is  a  solution  of  (2.8b)  defined 
in  R,  possessing  branchpoints  at  the  same  points  and  of  the  same  orders  as  g. 

i/'=  Im  [Pi(/)]  satisfies  equation  (2.8a)  and  can  be  interpreted  as  a 
stream  function  (in  the  A,  0-plane)  of  a  ( possible )  flow  pattern  of  a  com¬ 
pressible  fluid. 

Making  an  obvious  modification,  we  can  in  a  similar  way  extend  the 
definition  of  the  operator  P*  so  that  it  can  be  applied  to  functions  of  one 
variable  A  -f-  0  and  A  —  6  respectively,  thus  generating  solutions  of  equation 
(2.  12).  The  expression 

(3.25)  Pi(/i(A  +  0))  +P1(/2(A  — 0)) 

where  ft  and  f2  are  two  linearly  independent  functions,  will  represent  a 
(possible)  streamf unction  of  a  supersonic  flow  pattern. 

4.  The  integral  operator  of  the  second  kind  in  the  case  of  the  sim¬ 
plified  compressibility  equation.  The  integral  operator  of  the  first  kind, 
convenient  though  it  is  for  many  purposes,  has  the  disadvantage  that  it  does 

12  At  poles  and  logarithmic  singularities,  operators  P(f)  do  not,  in  general,  preserve 
certain  properties  of  P0{f)  which  are  essential  in  aerodynamical  applications;  in  these 
cases,  we  have  to  use  other  means  {see  [3,  9]),  in  order  to  produce  the  necessary 
singularities  of  4>. 
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not  represent  solutions  of  the  compressibility  equation  in  the  neighborhood 
of  the  sonic  line.  Furthermore,  it  has  the  disadvantage  from  the  practical 
point  of  view  that  the  PM  ( Z' ,  Z')  are  functions  of  two  variables  which 
makes  tabulation  of  the  values  of  P{n)  very  time-consuming. 

These  two  disadvantages  can  be  removed  by  the  use  of  another  operator — 
to  be  termed  “  operator  of  the  second  kind 99  ■ — for  which  the  P{n)  are  functions 
of  one  variable  only  and  which  yields  a  representation  of  the  streamfunction 
in  the  neighborhood  of  the  sonic  line.  This  operator  has  a  number  of  other 
distinctive  features  which  will  best  be  elucidated  by  the  detailed  discussion 
of  the  so-called  “  simplified  ”  compressibility  equation,  i.  e.,  where  Ar  =  Af 
=  (ISA)"1  in  equation  (2.8a)  or  alternatively,  where  F  =  Af  =  5/144A2  in 
equation  (2.14). 

According  to  Theorem  3.  1,  any  function  E  of  the  form 
(4.1)  E  =  HE*, 

where  E*  is  a  solution  of  the  equation 

(4.  2)  G2(E*)  =  (1  —  P)Ezt  —  (1  /t)E*i  +  %ZtE*Ei  +  2 ZtFE*  =  0, 
and  H  is  defined  by 

N(r)dT] 

■oo 

=  ( 1  —  M2)  [  1  +  i  ( h  —  1 )  M- 2  ]  -P2  * 

—  iST0( —  2A)-1'>((—  2A)2/3) 

with 

p((—  2A)2/3)  —  1  +  iS1( —  2A)2/3  +  ft(-~-  2A)4/«  +  .  • 

S0  ^  2(2fc+1)/<6fc-6)3“1/6(ik  +  1) 

(4.  4)  ft  =  (1/10)  (3/4 )2P(h  +  1)“V3  (2 Jc  +  5), 

ft - (1/1400)  (3/4) 4/3 (*^7  +  l)“2/3(64fc2  +  70&  +  75),  •  •  • 

may  be  used  as  a  generating  function  of  our  operator.  In  the  case  of 
the  simplified  equation,  we  have  Af  =  1/(12A),  p(( — 2A)2/3)  =1,  Hf(2A) 
=  ft  (—  2A)  -1/6,  Af  =  5/144A2. 

We  now  introduce  a  new  variable 

(4.  5)  u  =  t2Z/ ( Z  -f-  Z) , 

and  we  shall  show  that  in  the  simplified  case  there  exist  solutions  of  (4.  2) 
which  are  functions  of  u  alone. 
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Lemma  4. 1. 

(4.6)  £*f(A,M)  =  AxF(l/6,  5/6,  1/2,  —  t2(X  +  id)/—  2\) 

+  B1{~t2{\  +  iO)/—2kY*¥{2/%,  4/3,  3/2,  —  t2  (X  +  iO) /— 2\) 

(where  F(a,  /?,  y,  X)  denotes  the  hypergeometric  function  and  Ax  and  B± 
are  arbitrary  constants)  is  the  most  general  solution  of  (4.  2)  which  is  a 
function  of  u  alone.'  Let  us  note  that  there  exist  other  solutions  of  (4.  2) 
which  are  functions  of  one  variable,  see  (4.  20).) 

Proof.  We  shall  show  that  (4.2)  can  be  reduced  to  an  ordinary 
differential  equation  whose  solution  is  (4.  6).  A  formal  computation  yields 

du/dt  =  2  u/t,  du/dZ  =  (t2u  —  u2)/t2Z,  du/dZ  =  —  u2/t2Z, 

E*]z  =  —  u2t-2Z~'E*\Vj,  E*\Zt  =  —  2 uH-sZ-'[uE*U  +  E* fM] 

E*\zz  =  —  uH-*Z~2  [  ( ut 2  —  u2)E*]uu  +(t2  —  2u)E*  f«], 

Ff  =  (5/3  6)u2t-*Z~2. 

Substituting  the  above  expressions  into  (4.  2)  we  obtain 

G2(E* f)  =  —  2uH-sZ[u(l'—u)E*juu  +  (i  —  2u)E*%—  (5/36)F*f]  =  0. 

The  equation 

(4.7)  w(l  —  u)E*\tltl+  (1  —  2  u)E*%  —  (5/36)F*f  =  0 

is  a  hypergeometric  equation  whose  general  solution  can  be  represented  in 
the  form 

(4.  8)  F*f  =  ^(1/6,  5/6, 1/2,  u)  +  F^F(2/3,  4/3,  3/2,  u)  |  u  |  <  1 

(4.  80  =  A2u-V«¥(  1/6,  2/3, 1/3, 1/u ) 

+  B2u~5/6¥(5/6,  4/3,  5/3, 1/u),  \  u  |  >  1. 

Eeplacing  u  by  the  right-hand  side  of  (4.  5),  we  arrive  at  (4.  6). 

Thus,  combining  (4. 1)  and  (4.  6),  we  obtain  for  the  generating  function 
in  the  simplified  case  13 

13  We  note  that  in  many  instances  we  may  omit  the  second  term  on  the  right-hand 
side  of  (4.9)  since 

J\f(2/3;  4/3,  3/2,  — 12(\  +  id)/—  2\)f  (iZ(\  —  tz)  )dt/(l—  t*)l  =  0 

t=- 1 

if  f  is  regular  at  Z  =  0. 
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(4.  9)  Ef(X,  6,  t )  =  A1S0(—  2A)-1/^(l/6,  5/6, 1/2,  —t2( A  +  i$)/—  2A) 

+  W0(“  2A)-2/3[—  *2(A  +  t0)]*F(2/3,4/3,  3/2,  —  tf2(A  +  iO)/—  2A), 

|— *2(A  +  t0)./— 2A|  <  1 

(4.9')  =^o[(~^)(A  +  ^)]-1/6r(l/6,2/3,l/3,  —  2A/—  *2(A  +  i0)] 

+  £2£0(—  2A)_2/3[( —  t2)  (A  +  t<9)  ]-5/6F(5/6,  4/3,  5/3,  —2 A/—  *2(A  +  iO) ), 

|  —  2A/— tf2(A  +  10)  |  <  1. 

In  Theorem  3. 1  we  proved  the  existence  of  a  generating  function  by 
means  of  which  we  can  obtain  solutions  of  the  compressibility  equation  in  a 
sufficiently  small  neighborhood  of  the  origin.  We  shall  now  show  that  in 
the  case  of  the  simplified  equation  and  the  operator  of  the  second  type  this 
result  in  the  small  can  be  replaced  by  a  result  in  the  large.  This  result 
enables  us,  from  the  behavior  of  the  associate  f,  to  make  conclusions  con¬ 
cerning  the  behavior  of  the  generated  solution  of  the  compressibility  equation 
in  its  entire  domain  of  definition.  An  exact  formulation  of  this  statement, 
at  first  for  the  subsonic  region,  is  given  in  the  following  theorem: 

Theorem  4.1.  Suppose  that  the  function 

(4. 10)  g(Z)  =  f  /[*Z(1-P)]<H/(1  — *■)» 

J  *=-i 

is  regular  in  a  region  B  (situated  on  a  Riemann  surface)  which  possesses  in 
its  interior  a  finite  number  of  branchpoints /4  each  of  finite  order .  Let 
further  the  projection  of  B  on  the  schlicht  A,  9-plane  lie  in  E[ —  oo  <  A  <  0]. 

The  function 

(4.11)  if  —  Im[Pf2(/)], 

pf.(/)- f  z=x+ie 

^  c 

(where  C  is  a  suitably  chosen  curve  in  the  complex  t-plane  connecting  the 
points  t  =  —  1  and  t  =  1)  is  a  solution  of  (2.8)  with  N  =  Wf  —  1/12A; 
this  solution  is  defined  in  B  and  possesses  branchpoints  at  the  same  points 
and  of  the  same  order  as  (4.10).  if  =  Im[Pf2(/)]  can  be  interpreted  as  a 
streamf unction  (in  the  A,  0-plane)  of  a  (possible)  flow  pattern  of  a  com¬ 
pressible  fluid  (for  the  simplified  compressibility  equation). 

14  We  assume  here  that  the  only  singularities  of  g  in  B  are  branchpoints.  In 
applying  the  integral  operator  method  in  the  ease  where  g  has  poles  or  logarithmic 
singularities,  certain  modifications,  indicated  in  [9,  p.  469,  footnote  14],  are  needed. 
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Proof.  In  order  to  prove  our  statement,  we  have  to  show  that  by  (4.  6) 
and  (4. 10)  and  slight  modifications  of  these  formulas,  E\2(\,  0,  t)  is  defined 
for  all  values  A  <  0,  *  co  <*0  <C  °o  and  for  values  t  belonging  to  a  suitable 
simple,  sufficiently  smooth  curve  in  the  complex  tf-plane  which  connects 
£  =  and  f  =  l.  Obviously,  this  curve  has  to  avoid  the  points  t  =  0, 
^  ~  —  (»A/(X  +  *0))J  as  these  would  give  rise  to  singularities  of  the  hyper- 
geometric  function.  On  the  other  hand,  any  such  curve  will  be  suitable  for 
our  purposes.  However,  with  a  view  to  the  subsequent  generalization  of  our 
procedure  to  the  “  exact  ”  case,  we  shall  use  two  special  paths  of  integration, 
G1  and  02}  the  former  apart  from  its  terminals  t  =  ±  1,  inside  E[|  t  |  <-l], 
and  the  latter  outside  E[|  t  j  1].  Ci  will  be  used  for  values  A,  0  satisfying 
|  (A  +  i6)/2 A  |  <  1  and  C2  for  the  case  |  (A  +  iO) /2A  |  >  1. 

It  should  be  noted  that  the  two  terminals  of  the  integration  path,  viz., 
t=±  1,  will  never  be  singularities  of  Ej2(A,0>t)  since,  for  t  =  hh  1, 
(A  -f-  i6)/2 A  1  for  real  A  and  0. 

The  expressions  thus  obtained  will  not  necessarily  be  analytical  con¬ 
tinuations  of  each  other  (qua  functions  of  A,  0).  Since,  however,  the  hyper¬ 
geometric  equation  has  only  two  linearly  independent  solutions,  the  constants, 
A i,  and  A2,  B2  can  always  be  so  adjusted  as  to  make  these  two  solutions 

analytical  continuations  of  each  other. 

Remark.  It  would,  of  course,  also  be  possible  to  characterize  the  path 
of  integration  in  a  manner  which  is  topologically  invariant  with  regard  to 
the  way  the  singular  points  of  Ef2(A,  0,  t)  are  by-passed;  using  this  definition 
we  would,  for  any  value  (A,  0),  obtain  one  and  the  same  function  ^(A,  0) . 
However,  although  this  procedure  has  some  theoretical  advantages,  its  actual 
carrying. out  may  give  rise  to  certain  difficulties;  in  practical  applications 
it  is  much  easier  to  assure  analytical  continuation  by  the  determination  of 
the  constants  A2)  B2  if  A1}  B±  are  given,  or  vice  versa. 

The  integral  representation  (4. 11)  can  be  immediately  generalized  to 
the  supersonic  case  where  it  will  produce,  in  an  analogous  manner,  solutions 
of  (2. 12)  with  N1~Nj[1  (see  (2. 18)).  Indeed,  replacing  A  by  the  variable 
w  =  A  +  tA  and  considering  the  solution  if/(<o,0)  of  (2.8a)  in  the  plane 
A  =  0,  it  is  seen  that  if/(iA,0)  satisfies  equation  (2.12)  with  N=  1/12A. 
Repeating  the  procedure  which  led  to  the  generating  function  (4.  9)  in  the 
subsonic  case,  we  now  obtain  the  generating  function 

(4. 12)  tff  =  (diS0/ (2A)1/6)F(l/6,  5/6, 1/2,  t2( A  +  0)/2A) 

+  (Mo(A  +  0)H/ (2A)2/3)F(2/3,  4/3,  3/2,  t2(A  +  0)/2A), 

|*2(A  +  0)/2A|<l. 
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If  F(<x1?  /31}  yi,X)  denotes  only  the  hypergeometric  series  and  not  the  hyper¬ 
geometric  function,  (4. 12)  has  to  be  replaced,  for  |  t2(A-\-  0)/2A  |  >1,  by 

(4.13)  Ef  =  (a2S0/(t2(A  +  0))1/6)T(l/6,  2/3, 1/3,  2A/*2(A  +  0)) 

+  (&2$0(2A)2/3/(tf2(A  +  0)  )5/6)F(5/6,  4/3,  5/3,  %A/t*(A  +  6)), 

where  the  constants  a2,  b2  are  easily  expressible  in  terms  of  ax,  bx. 


If  e  =  A  and  t2  =  1,  there  arise  certain  difficulties,  since  the  hyper¬ 
geometric  functions  in  (4. 12)  will  then  become  singular.  By  the  trans¬ 
formation  formulas  of  the  hypergeometric  function,  (4. 12)  may  be  written 
'in  the  neighborhood  of  ( (A  +  6)/2A)t2  =  1,  in  the  form 

(4.14)  Ef  =  (^o/(^A)1/6)F(1/6,5/6,3/2,  1  -  *2(A  +  0)/2A) 
b3S6 (2A/  (2A  -  t2  ( A  +  6)  ) )  ^  (1/3,  —  1/3, 1/2, 1  -  t2  ( A  +  0)/2 A) . 

In  order  to  avoid  the  complications  which  arise  from  the  fact  that  the  second 
term  of  (4. 14)  has  a  singularity  for  A  =  0,  t2  =  1,  we  shall  therefore  take 
b3  =  0.  The  function  Ef  will  accordingly  be  of  the  form 

(4.15)  A'f==(^0/(2A)1/6)F(l/6,5/6,3/2,  1  —  tf2(A  +  0)/2A). 

We  note  further  that  all  these  considerations  can  be  repeated  with  A 
replaced  by  — A.  Our  operator  will  therefore  yield  two  independent  types 
of  solutions,  depending  on  whether  the  argument  of  the  associate  function 
is  taken  as  A  -f-  6  or  A  —  0. 

The  exact  conditions  under  which  our  operator  can  generate  solutions 
of  the  compressibility  equation  in  the  supersonic  case  -  are  given  in  the 
following  theorem : 


Theorem  4.2.  Suppose  fs(£),  *  — 1,2,  are  real  functions  of  the  real 
variable  {  and  everywhere  differentiable  with  the  possible  exception  of  C  —  0 ; 
suppose  further  that  in  a  fixed  neighborhood  of  £  =  0,  fs  can  be  approximated 
to  any  prescribed  degree  of  accuracy  by  the  expressions  of  the  form 

!,>  1.  Then 

(4.16)  ^(A,  6)  =  Rfi(fi)  +  R\z(f2)} 


represents  a  solution  of  the  compressibility  equation ,  which  is  defined  for  any 
A  >  0  and  can  be  interpreted  as  a  streamfunction  of  a  (possible)  supersonic 
flow  pattern. 
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Proof.  By  (4. 12)  we  have 
(4.17)  Rf1(/1)=«i/(2  A)V» 

X  5/6’  1/2>  <2((A  +  0/8A )/(  (A  +  9)  (1  —  P) /2)  (1  —  t‘)-*dt 

+  c1(A  +  0)i(2A)-2/s 

X  JjF(3/3, 4/8,  8/2,  <2( A  +  0)/2A )/( (A  +  9)  (1  —  i2)/2)  (1  —  <2)-i^. 

In  view  of  A  >  0,  the  only  values  of  (A,  6)  for  which  this  expression  may  not  be 
differentiable  are  those  for  which  A  -f-  6  —  0  or  for  which  u  =  t2(A.  -{-  0)  /2A 
coincides  with  either  of  the  values  0, 1,  oo  for  —  1  <  t  <  1.  The  case  w  =  oo 
is  ruled  out  because  of  A  >  0 ;  the  case  u  =  1,  although  corresponding  to  a 
singularity  of  the  hypergeometric  equation,  does  not  give  rise  to  a  singularity 
of  R ,  by  virtue  of  our  particular  choice  (4. 15)  of  Ef.  The  case  u  =  0  need 
only  to  be  considered  for  A  -}-  0  =  0,  as  t  —  0  obviously  does  not  cause  any 
difficulties. 

Under  our  assumptions,  it  is  sufficient  to  consider  the  special  case 

(4- 18)  fi(t)  *  >  1, 

in  which  (4. 17)  reduces  to 

(4.19)  RUiifi,*)  =a1/(2Ay/Q 

X  Jf(  1/6,  5/6, 1/2,  A  +  0)/2A)  (A  +  0Y(t  —  t*y-idt 

+  Ci(2A)-2/3  J  <F(2/3, 4/3,  3/2,  <2(A  +  0)/2A)  (A  +  0)^(1  —  t^dt. 

For  A  -j-  6  0,  this  expression  is  obviously  continuous.  The  same  is  also 

true  of  the  derivative  dRj2tl/dJif.  Indeed,  0Rf2il/mf  (see  (2.18))  behaves 
in  the  neighborhood  of  A  -f-  0  ==  0  like  (A  -j-  S)^1 ;  in  view  of  k  =>  1,  it  there¬ 
fore  remains  continuous  there.  This  completes  the  proof  of  Theorem  4. 1. 

Before  we  proceed  to  investigate  the  behavior  of  our  solutions  on  the 
sonic  line,  we  note  that  (4.  6)  is  not  the  only  solution  which  depends  only 
on  one  variable.  If  we  write 

(4- g0)  p  =  t2(X  —  A0  -f-  i8) /2\ 

and  try  to  solve  (4.  2)  by  a  function  of  one  variable  /*,  a  formal  computation 
shows  that  E"\  as  a  function  of  fx  must  satisfy  the  hypergeometric  equation 
(4.  7)  with  u  replaced  by  fi. 
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Accordingly,  the  generating  function  (4.  6)  for  the  operator  of  the  second 
kind  may  be  replaced  by 

(4.21)  A'f  (A,  0,  t)  =  A±( — 2A)1/6F(l/6,  5/6, 1/2,  *2(A-A0  +  i0)/2\). 

Besides  its  greater  generality,  the  generating  function  (4.  21  with  A0^0 
has  a  number  of  additional  features,  which  make  it  superior  to  (4.  6)  in 
many  cases.  In  (4.6)  the  point  A  =  0,  0  =  0  is  a  singularity  since,  by 
approaching  this  point  in  a  suitable  manner,  the  argument  of  the  hyper¬ 
geometric  function  can  be  given  an  arbitrary  value.  In  (4.21)  such  a 
singularity  does  not  exist,  since  A  and  A  —  Ao  4-  iS  cannot  vanish  simul¬ 
taneously  if  A0  =^=  0. 

Another  singularity,  which  can  be  removed  by  using  (4.  21)  with  A0  0 
instead  of  (4.6),  occurs  in  the  supersonic  case.  For  t2  =  1  and  A  =  0, 
the  second  term  of  the  right-hand  side  of  (4.  14)  becomes  singular.  In  order 
to  allow  for  this  case,  we  had  to  assume  that  lz  =  0,  thus  somewhat  restricting 
the  generalitv  of  the  solutions  we  could  obtain.  Setting  A  =  'iA,  it  is  seen 
that  in  the  supersonic  case  the  argument  of  the  hypergeometric  function  in 

(4.21)  becomes 

(4.22)  ju,  =  £2(A  +  i\0  +  0)/2A. 

For  t  =  db  1,  we  have  /*,  ^  1  for  real  values  of  A,  0;  the  singularity  in  question 
is  therefore  removed  and  the  constant  &3  in  (4.14)  may  now  take  any 
arbitrary  value.  Moreover,  (4.22)  shows  that  ft  =^=0  for  A-j-0  =  O;  as  a 
consequence,  the  line  A  -f-  0  =  0  loses  its  singular  character  and  the  dis¬ 
cussion  of  the  equivalent  of  (4. 17)  is  considerably  simplified. 

Thus,  by  the  use  of  integral  operators  of  the  second  kind  with  the 
generating  functions  (4.  9)  and  (4. 12)  as  well  as  (4.  21)  we  obtain  solutions 
^Im [p2(g)]  which  are  defined  in  the  subsonic  and  supersonic  regions 
respectively.  If  g  is  defined  in  a  (not  necessarily  schlicht)  domain  G\  in  the 
subsonic  region  and  has,  as  its  only  singularities  in  G1}  branchpoints  of  finite 
order  (but  not  poles  or  logarithmic  singularities),  then  the  generated  function 
is  again  defined  in  Gh.  (In  particular,  it  has  branchpoints  at  the  same  points 
and  of  the  same  order  as  g). 

If  g  is  a  twice  differentiable  function  of  one  real  variable  in  a  schlicht 
domain  G2  in  the  supersonic  region,  the  generated  function  will  be  a  solution 
of  (2.18). 

Thus  by  this  procedure,  we  can  generate  solutions  of  (2.  2)  which  are 
defined  in  adjacent  domains,  one  in  the  subsonic,  the  other  in  the  supersonic 
region.  Our  object  is  to  show  that  these  solutions  can  be  continued  across 
the  sonic  line. 
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This  fact  is  immediately  seen  if  we  introduce  a  new  variable/5 

(4.23)  S  =  ^  for  A  <  0, 

s  =  —  A2/3  for  A  >  0. 

The  variable  s  =  s(i!f|)  considered  as  a  function  of  Mach  number  Iff, 
possesses  the  property  that  s(l)  =0,  and  that 

(4.24)  ds(Mj)/dM\  = —  25/3(3&  +  3)_2/3Mf  +0(1  — Iff2) 

is  non-vanishing  and  bounded  in  a  sufficiently  small  neighborhood  of  Aff  =  1. 

If  A0  +=  0,  and  if  g  is  regular  for  (—  60  <  0  <  0o,  s  =  A  =  A  =  0)  then 
the  generated  function  is  an  analytic  function  of  0  and  5  (and  therefore  of  0 
and  Aff)  for  (  0o  If:  0  <  0O}  — 5(0)  <1  s  <  s(0)),  s<°)  sufficiently  small. 

If  A0  =  0,  we  have  to  assume  that  lim  Z~5/6f(Z)  exists  (or  alternately 

that  Hm  Z1!Gg  (Z)  exists)  in  order  to  assume  that  the  generated  function 
\f/(0,s)  is  regular  also  at  the  point  6  =  0,  5  =  0. 


5.  Integral  operator  of  the  second  kind  in  the  case  of  the  “  exact  ” 
compressibility  equation.  As  we  mentioned  before,  any  solution  E*  of  the 
equation 

(5. 1)  JE*£t  —  t2E*zt  —  (1  /t)E*i  +  2  ZtE*zz  +  2ZtFE*  =  0, 


XZ+2 

Ndr],  (see  (4.3)),  is  a  generating  function  of 

00' 

an  integral  operator  P (/)  (see  (3.4b)  which  produces  solutions  of  equation 
(2.8b).  The  series  1  +  2  (t2Z)nQ(-n)  (2A),  where  the  (+n)?s  are  solutions 

n= 1 

of  the  system: 


s%  2X 

(5.2)  +^x(«)+4Fe<»)=0,  — 4  f  Fi 

*S  -OO 


(see  (87)  of  [6]),  is  a  solution  of  (5.1).  The  above  series  converges  for 
|  ^  |  <  2  |A  |  and  reduces  to  the  first  summand  of  (4.  6)  for  F  =  F\. 

It  is  therefore  desirable  to  obtain  solutions  of  (5. 1)  which  are  defined 
in  |  Z  |  >  2  |  A  j  .  In  this  section  we  shall  determine  two  power  series 
(^  —  2,  2),  both  converging  in  j  Z  j  >  2  |  A  |  ,  such  that  in  the  simplified 
case,  A2EE*(1)  +  B2RE*(2)  reduces  to  (4.9'). 


15  We  note  that  in  the  simplified  case  we  have 

S  = - 2  ■  3-2/3  {  2/ {&  —  1)  )  (2-*)/(3*-3)Hf. 
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Theorem  5.1.  Let 

(5.3)  g(»,.)=|o»(”''>(-A)^*Wsl(w)  « =  1, 2, 

P=0 

a  of  functions  which  are  connected  by  the  relations 
(5.  4a)  q\\{0’K)  +  —  0 

(5.  4b)  2(ft  +  }*)gx(,wr)  +  gx\(w+1,'f)  +  4:Fq(n+1'K)  —  0, 

71  — —  1,  2,  ■  *  *  j  K  lj 

Theft  each  of  the  functions 

00 

(5.5)  =  2  g(n,K)/ (—  t2Z)n~*+w3)K 

is  a  solution  of  (5. 1).  Each  of  the  series  converges  in  E[ 2  \  X  |  <  |  Z  |  ]. 

Proof.  Substituting  the  series  (5. 5)  into  (5. 1)  and  equating  the 
coefficients  of  H6^>/3(— Z)-<6«+1>/6,  ft  =  —  1,  0, 1,  •  •  • ,  to  0,  we  obtain  the 
following  set  of  equations: 

(5.  6)  qzz{0,K)  +  Fq(°’K)  =  0,  (ft  +  §*)gz (n,/c)  +  qzz[n+1,K)  +  Fq{n+1’K)  =  0, 

ft  =  0, 1,  •  *  • , 

which,  if  we  assume  that  the  q^K)  are  functions  of  A  alone,  result  in  the 

equations  (5.4a)  and  (5.4b). 

According  to  {2. 17),  F  can  be  written  in  the  form 

(5.7)  F  =  s_3#(s),  s=( — A)2/3, 

where  S(s)  =  a0  +  +  a2s2  +  •  •  * ,  «o  =  5/144,  c?i  =  0,  •  •  • ,  see  (2. 17), 

considered  as  a  function  of  5,  is  regular  in  a  circle  of  radius  s0  say,  with 

center  at  the  origin. 

Introducing  the  variable  5,  we  obtain 

gx(«> _ (ds/d(—\))(dqM/ds)  — *(— A  )-'/*(dq^/ds) 

=  —  s~%(dq  (n’K)/ds) 

gxx(n+i,«)  _  (4/9)  [—  fr-% +  r1?,/"*1'’1]. 

Hence,  the  system  (5.4a),  (5.4b)  assumes  the  form 

(5.  9a)  s2qss (0,l!>  —  isqs(0^  +  9S(s)q<°’K>  =  0 

(5.9b)  _3(«+Sk)s5/22/»'‘>  +  —  Isq*^1’*  +  9S(s)q^1'K)  =  0; 

n  —  0, 1,  2,  3,  •  ■  • . 
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Lemma  5. 1.  There  exists  a  system  of  solutions  q <»•*>  n  —  0  1  2  ■  ■  ■ 
of  (5.  9a),  (5.  9b)  of  the  form 

(5. 10)  jM_s(WW(Wl)„riW  {s)>  TkM{s)  =  |  Co<M)  _^0 

J/=0 

where  each  2V«>  (a)  is  a  power  series  which  converges  in  the  circle  |  s  |  <  *0. 

Pr°°f.  We  shall  first  prove  the  above  lemma  for  q  (°*>  and  then  for 
arbitrary  n  by  induction  on  n,  i.  e.,  we  shall  show  that  if  the  lemma  holds  for 

then  dt  must  for  Let  us  first  consider  the  homogeneous 

equation 

•  (5-n)  S2wss  —  ±SWS  +  9S(s)w  =  0, 

The  indicial  equation  (see  [15],  p.  225)  is 

(5.  12)  p(p  1)  —  ip  5/16  =  o,  or  Pl  =  5/4,  p2  =  1/4. 

By  substituting  w  —  sx^y,  we  obtain  the  equation 

(5-  13)  y"+  9s-*[S(s)  —  5/144]*,  =  0. 

Since  5  ~[S(s)  5/144]  is  a  regular  function  of  s  for  |  s  |  <  s0,  (see  5.  7)) 

we  may  choose  the  two  particular  solutions  of  (5. 13)  as 

=  Cr0(0’1)  +  O2(0>1)52  +  C3(0'1)s3  +  *  •  • 
y*  =  C0i0’2)s  +  C^O’^s2  -f  O2(0’2)53  +  •  •  . 

which  yield 

.  -==sV*wM(s) 

_  (70«>.2)SV4+  0^0,21^9/4  _|_  C2(04)s18/4  +  .  .  .==^^(2)  ^ 

In  the  case  of  the  simplified  equation  we  have  [see  (4.  9'),  (4. 1),  (4. 3)] 
?t(0,1)  =  CV0’1)  ( — A)1''6,  Cfai”'1)  =  2'/°, 

gt< M)  =  Cfo*0'2*  (—  A)5/«,  C*f  <».2>  _  2»/«. 

In  order  to  obtain  analogous  series  for  the  “  exact  ”  ease  we  choose 
(5. 14)  Co*0'1)  =  Cf/0'1’  =  21/6,  (70(o-2)  =  Cf„<0’2>  =  a»/». 

Let  us  now  instead  of  (5.9a)  consider  the  solution  of  the  non-homo- 
geneous  equation  (5.  9b)  and  let  us  assume  that  for  n  >  0,  we  have  already 
proved  that  has  the  form  (5.  10)  and  !,<•)(«)  converges  for  |  s  |  <  So. 

In  order  to  prove  that  there  exists  a  solution  of  the  equation  (5  9b) 
of  the  form  _,«/»<•*«/»>  .>r,<»«)  («),  we  proceed  as  follows- 
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Let 
(5. 15) 


q(n+i,ic)  =  w(s)'lt($), 


where  w(s)  =  sK-s/4TT(/c)  (s)  is  a  solution  (see  above)  of  the  homogeneous 
equation  (5.11).  Then  u  will  satisfy  the  equation 

(5.16)  •  wusa  +  (2ws  —  w/2s)us  =  3(n  +  ^k)s%^k). 

The  particular  solution  of  this  equation  is  easily  verified  to  be 

(5.17)  m  =  3(«  +  $k)  f°w-2st(j‘°wqsMds)ds 

Expanding  qln*1’,:> ,  we  obtain 

(5.19)  g  <*♦»•*>(«)  '(*), 

2V”+1>  (0) 

=  [(6»  —  3  +  4k)  (3ft  +  2k)/(3»  +  3)  (3n  +  4k  —  3)]ZV(’*)  (0)  ^  0 
the  series  for  2V’*H'1)  (s)  converging  for  |  s  |  <  -V 

This  completes  the  proof  of  Lemma  5. 1.  In  order  to  prove  Theorem  5. 1, 
it  remains  to  be  shown  that  it  is  legitimate  to  interchange  the  order  of 
summation  in  the  double  series 

(5.  20)  2  (—  t2Z)-(n~*+{2/s)K)  2  Gvin,K)  ( — A)b^+(2/3)  {k+v) 

'  n=o  '  *=0 

for  2  |  A  |  <  |  Z  |  .  For  this  purpose,  we  shall  prove  10 


Lemma  5.  2. 
(5.  21) 


(7Jbt{n+1’1)  |  g  2W+1M/V 


for  where  C^w+1,1)  are  the  coefficients  of  the  series  (5.2). 

Here,  p0  and  M  are  sufficiently  large  constants,  and  Si  =  s0(l  +  0-1- 

Proof.  We  shall  give  a  proof  by  induction.  Consider  at  first  the  CVn,1)> 
n  __  2,  •  •  *.  If  we  substitute  the  power  series  (5.3)  into  (5.4b),  it 

becomes  evident  that  the  C0(B,1)  depend  only  upon  <x0,  and  are  independent 
of  the  remaining  coefficients  a»,  n  >  0,  of  the  series  expansion  of  8(s), 
(see  (5.7)).  On  the  other  hand,  if  we  substitute  a0  =  5/144,  <xn  =  0  for 
n>  0,  we  obtain  the  simplified  case  considered  in  (4.9'),  where  we  had  a 
representation  for  E\  as  a  power  series  in  (2A/ — t2Z).  Using  the  fact  that 

16  We  shall  consider  here  only  the  case  when  k  =  1.  Exactly  the  same  proof  holds 
for  k  =  2. 
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2£*f  =  Ej/Hj  =  $0_1( —  2A)1/6i?t  (see  (4. 1)  if.)  we  obtain  two  power  series 
for  jE^f ;  the  first  of  which  corresponds  to  the  case  k  =  1  and  the  second  to 
the  case  k  =  2.  Thus,  setting  A2  =  1  in  (4.  9'),  we  have 

(5.22)  (—  2A/—  ^Z)1/^  (1/6,  2/3,  1/3,  —  2A/—  *2Z) 

=  (— 2A/—  t2Z)1/Q 

,  £&(&  +  !)  ■  '  ■  (i  +  w)  "i(S  +  1)  (f  +  n)  (  o\  /  >27\„+l/6 

whence 

(5.  33)  Co*0’1*  =  Cfo'0'1’  —  21/e,  00<n’1)  =  CV’1* 

Ki  +  1)'  -  •d+ft)(f)(f  +  1)-  ’  ■(■!  +  n)2»»V 

J(s  +  1)  ‘  *  '  (i  +  n)  * n- 

=  1,  2,  3,  *  *  * , 

from  which  the  inequality  (5.21)  for  /x=0  follows. 

Remark.  Since  in  this  case  y  =  /x  =  0,  the  number  s0  in  (5.  21)  where 

v  is  replaced  by  /x)  may  be  given  any  positive  value.  For  n  =  —  1,  and  an  arbi- 

00 

trary  /x,  the  inequality  (5.  21)  follows  from  the  fact  that  g(0,1>  =  2  0  p.i0,1)  s^4 

00  PL-  o 

is  a  solution  of  (5.  9a)  and  the  fact  that  S(s)  =  2  V*  is  regular  in  the 

fi=Q 

circles  of  the  radius  s0  (see  also  (5. 13)) ;  accordingly, 


(5.  24) 


holds,  if  T  is  a  sufficiently  large  constant. 

We  now  proceed  to  the  proof  (by  induction)  of  (5.21)  for  — 1 
and  /x  >  0.  Let  us  assume  that  this  inequality  holds  for  some  n  +  1,  and 
fj,  <  v  —  1  as  well  as  for  N  <  n  and  /x  ^  v  +  1.  We  shall  prove  that  (5.  21) 
then  holds  for  N  —  n  +  1,  /x  =  v. 

If  we  substitute  the  series  (5.3)  (with  k  =  1)  into  (5.4b),  we  obtain 
the  relation 

—  »(»  +  $)(»  +  *  +  b)Ov^  +  (»  +  §  +  fv)  (n  +  i  +  f v)<V«w> 

(5.  25)  v-i 

V  '  +  4a0(7v(”+1,1)  +  2  4(y*+1’1)aV-M  =  0. 

/*= 0 

Hence, 

(5.26)  |  Cfv<"tl'»[l+  (4.«o/(»+5  +  §>')(«  + J'  +  fv))]  I 

5S  2»+1M/Vi’"1{[(1  +  2/3«)/ (1  +  V6»  +  (S/3)  (v/») )] 

+  4vr/(«  +  f  +  fv)  (n  +  J  +  Jr) }. 
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If  n  +  fv  >  p0,  then 

(5.  27)  |  Cv <n+1’V  [1  +  4a0/p02]  |  g  2n+1M/s0s1v~1[l  +  6T/p0], 

or 

|  Cv |  ^  2n+1M(l  +  c)/viM  =  2”+1M/$/. 

This  completes  the  proof  of  Lemma  5.  2. 

In  order  to  show  that  it  is  legitimate  to  interchange  the  order  of  summa¬ 
tion  in  (5.  20),  we  shall  show  that  for 

|  A  | 2/3  g  $o  —  €;  |  —  2A/ —  t2Z  |  <  1  —  c,  where  €  >  0, 

the  series  converges  absolutely  and  uniformly.  Indeed,  by  Lemma  5.  2  we 
obtain 

00  oo 

(5.28)  |2  ( t2Z )  -W*)  -»  2  v (n>1)  (—  A)  n+{1/6> +  (2/3>v  | 

«,=0  j/=0 

00  00 

g  M  2  2  2^— A)”+(1/6)  +  (2/3>VV  |  tf2Z  |  <V6)+n 

n=0 

=  M  |  —  2 A./ —  t2Z  |  V«[i  _  |  __  2A/—  t2Z  |]-i[l  —  |  A2/y5l|  ]-i 

which  for  |  A  |  2/3  <  s0  —  e,  |  —  2 A/ —  tf2Z  |  <  1  —  c,  0  <  c  <  1,  becomes 
smaller  than  M(1  —  c)1/6/e2,  which  shows  that  the  series  converges  absolutely 
and  uniformly. 

In  order  to  obtain  a  continuation  of  a  given  streamfunction  to  the 
supersonic  region,  i.  e.,  for  imaginary  values  of  A,  it  is  convenient  to  replace 
A  by  the  variable  s  —  ( —  A) 2/3  used  before  in  a  different  context  (see  (4.  23) ). 
Then  the .  generating  function  in  the  subsonic  case  may  be  written  as 

(5.  29)  J?<*>  (A,  0,  t)  =  H(2X)E*W  (A,  0,  t) 

OO  00 

—  802-1/ep(s)sK-1 2  2  Cv^sW^Vs'’ / (— Z)«-i*<2/3)« 

«=0  y=o 

|  5  I  <  So,  8  I  *  |  3/2  <  I  z  | ,  *  =  1,  2, 

(see  (4.3)  and  (5.5))  where  now 

(5.29a)  Z  =  —  s3/2  +  ^0. 

By  replacing  A  by  iA,  we  see  that  ( — A) 2/3  is  changed  to  • — A2/8.  Thus 
we  obtain  as  the  generating  function  in  the  supersonic  case 

(5.30)  EW(A,0,t) 

00  oo 

=  %-1^8<,p(s)sK-1 2  2  inCv(n’^  (—s)  u 

|  s  |  <  so,  3  |  *  |  3/2  <  |  Z  |  , 

where 
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(5.30a)  Z  =  i(—  s)V2  +  i6. 

Let  us  now  restrict  ourselves  to  a  neighborhood  of  the  point  A  =  0,  0  =  60 
(60 0)  lying  entirely  in  the  domain  D  =  E[|  s  |  <  s0,  \  s  |  <  3-1/3 1  q  |  2/3] . 

For  ( s ,  6)  e  E[s2  +  (0 —  0o)2  <  0o2/25],  we  have 
(5.31a)  (^0)t[1_(s3/2  —  i(O  —  0o))/iOo]v 

=  2^(0 - 60)mS{Z/2)n  for  M  <  1, 

m,n 

(5.31b)  Zv=  (i0„)'v[l  +  ((— s)3/2  +  (0  — 0„))/0]? 

=  2  &m»(0  —  00 ) "‘s' 3/5,11  for  M  >  1 

m,n 

because  then  both  |(s3/2' — i(9  —  0o))A0o|  and  |(( — s)3/2-|-i(0 —  0O))/0O| 
are  less  than 

(Ms/2  + 1 0— 0„  |)/|  | 

<  (3-1/2  |  0  |  +  |  0  — 0„  \)/\  0O  |  <  (2-3^  +  l)/5  <  1. 

Here-<imji  and  bmn  are  suitably  chosen  constants.  Thus  both  the  subsonic 
generating  function  and  the  associate  function,  and  therefore  P[/( — s3/2 
+  i0)],  may  be  expanded  in  integral  powers  of  (0  —  0O)  and  s1/2.  Similarly, 
in  the  supersonic  case,  P[/(i(—  s)3/2  + 10)  ]  may  be  expanded  in  integral 
powers  of  (0  —  0O)  and  ( — s)1/2.  Both  functions  are  determined  and  equal 
to  each  other  for  s  =  0. 

We  shall  show  that  if  f  is  regular  in  D,  the  coefficients  of  sv/2,  v  odd, 
vanish,  and  therefore  P[/( — s3/2  +  i0)]  is  an  analytic  function  of  s  and 
(0  —  0o)  at  (O,0o). 

Theorem  5.2.  If  f  is  regular  in  the  domain  D  =E[|s|  <  s0, 
|  s  |  <  3_1/3  |  6  |  2/3],  then  the  solution  P[/( — s3/2  -j-  iO )]  is  a  regular  function 
of  s  and  (6  —  0O)  in  D. 

Proof.  We  have  shown  that  P(/)  may  be  expresed  as  a  power  series 
in  s1/2  of  the  form 

00 

(5.32)  p(/)  -2M*  —  0o)W2. 

n=  0 

oo 

Using  the  variable  s,  and  noting  the  fact  that  N  =  (l/12)s~s/2(l  +  2  flWSw+1) 

n~  0 

(see  (2.16))  equation  (2.8a)  becomes 

00 

(^*  8  )  4i p8S  -j-  9s\p6d  -J-  2\j/s  2  Q'n+ i5n  =  0. 

n=o 

By  substituting  (5.  32)  in  (2.  8')  we  obtain  the  recursion  formulae: 
(5.33a)  ft1  =  fc3  =  0 
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(5. 33b)  hn+4 

(n- D/2 

==' —  [1/ (w  -f-  4)  (%  2) ]  [9^w_2  -f-  2  2  ( (3/^)  +  /)a(»/2)+i-i^2j+3] 

3-0 

for  n  odd; 


(5. 33c)  hn« 


n/2 


-  [1/  (fl  -f-  4)  (?&  -f-  2)]  [9h"n-Z  -j-  2  2  (1  y )  W/2)+l— J^2  j+2] 

i=o 

for  w  even. 


Since  hn,  where  n-  is  odd,  depends  only  on  the  previous  hm  for  odd  m, 
we  see  from  (5.  33a)  that  hn  —  0  for  all  odd  n. 

Thus  (5.32)  becomes 

00 

(5.  32')  P[/(— s8/2  +  t0)]  —  —  «o)s". 

Since  5  is  unchanged  by  the  result  of  putting  iA  for  A;  we  obtain 

00 

(5.  34)  P[/(f  (— s)3/2  +  iO)]  =  2  fn(0  —  0o)s\  . 

n- 0 

The  expressions  P[/( — s3/2-f-^)]  and  P[jf(i( — 5)3/2  +  ^)L  qua  functions 
of  5  and  (0  —  0O),  are  analytic  continuations  of  each  other  across  the  sonic  line. 

Thus,  assuming  that  the  associate  function  is  regular  in  a  sufficiently 
large  domain,  and  applying  the  integral  operator  of  the  second  kind,  we 
obtain  solutions  of  a  compressibility  equation  defined  in  four  adjacent  domains, 

Dx  —  E[Jf  <  1,  0  >  31/2  |  \  (M)  |]  +  E[If  >  1,0>  A(jJf)], 

'  D2  =  E[if  <  1,  |  0  |  <  31/2  |  \(M)  |], 

Z>3  =  E[1T  <  1,0  <  —  31/2 1 \(M)\]  +  E [M  >  1,0  <  —  3A(iH)], 

D4  =  B[M  >  1,  — 3A (M)  <  0<  A (M)]. 

(See  fig.  1,  p.  859.)  The  solutions  defined  in  Dx  and  D8  were  derived  in  the 
present  paper,  while  those  defined  in  D2  and  Z>4  were  derived  in  [6,  §  11]. 

In  the  simplified  case,  using  the  theory  of  hypergeometric  equations,  it 
was  possible  to  combine  these  representations  into  one,  which  yields  solutions 
of  (2.  2)  defined  in  the  whole  M-0- plane.  In  the  exact  compressibility  equation 
the  problem  remains  of  combining  these  four  representations  into  one.  This 
problem  can  be  attacked  by  using  the  integral  operator  of  the  first  kind  in 
addition  to  that  of  the  second  kind,  and,  in  analogy  to  the  simplified  case, 
developing  a  theory  of  differential  equations  with  singular  coefficients,  which 
would  furnish  us  with  information  corresponding  to  that  used  in  the  simplified 
case.  As  will  be  shown  elsewhere,  the  methods  of  the  Euchs  theory  for  ordinary 
differential  equations  can  be  generalized  to  the  case  of  partial  differential 
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equations  of  type  (2. 14)  with  F  being  an  analytic  function  of  two  complex 
variables  and  possessing  certain  singularities.  In  particular,  if  the  singularity 
surfaces  are  linear,  the  following  results  which  will  be  proved  in  a  subsequent 
paper,  will  be  valid: 

Theorem  5.  3.  Let  the  coefficient  F  of  $  in  equation 
(5.  35)  ( d2if//dZ1 2)  -j-  (d2ij//dZ22')  -p  Fip  —  0,  Z1  =  A.  -{-  tA,  Z 2  =  6  -f*  i0? 

have  the  form 

00  00 

(5.  36)  =  H  Am-1>n-iW’n-'Zn-\  =£  0, 

m- 0  n- 0 

W  —  Z1+iZ2,  Z  =  Z1  —  iZ2 
where  the  coefficients  A.m>n  satisfy  the  inequalities 
(5.  37)  |  Am-i,n-i  |  A/pinpon,  A  ^  1,  pi  0,  pz^>  0, 

being  suitably  chosen  constants. 

If  r  0  is  a  complex  constant  ivhich  satisfies  the  inequalities 

(5.38a)  r=^= — (m/2)  [1  dr  (1  + 

0  <  arg  (1  +  A.^/mn)*  <  ir/2, 

(5.38b) 

(5.38c)  —  tn 

for  m  =  0, 1,  •  •  •  ;  n  =  0, 1,  •  •  •  ;  (m,  n)  =^-(0,  0),  then  the  expression 

00 

(5.39)  1 f  = 

n- 0 

will  represent  a  solution  of  (5.35).  Here  B00  is  an  arbitrary  complex 
constant ,  and  the  Bmn  have  to  be  determined  from  the  equations 

(5.  40)  4[mn  +  rn —  ( J._lj_1w/4r) ~\Bmn 

m- 1  n  «- 1 

= - 2  2  Am_ju-i,ft-v-i-Sjuv - 2  A_itn-v-il5)nv* 

*>=0  V-Q 

(5.  38)  converges  in 

E [ |  W  |  <  pi/A.B,  |  Z  |  <  po/AB]  where  B  =  max  [1,  |  500  |',  M/4], 

M  =  max  { (mn  -j-  m  +  n)/|  mn  +  rn  —  4_1r_1mJ4_if_i  |}, 

m  =  0, 1,  2,  •  •  • ,  n  =  0, 1,  2,  •  ■  * ,  (mn)  =^=  (0,  0). 


and 
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Theorem  5.4.  Let  the  coefficient  F  of  if/  in  (5.35)  have  the  form 

00 

(5.41)  F  —  Zx~2A-2,q  7i{~xA- i(Z2)  "h  2  -4»(Z2)Z1n 

♦t= 0 

where 

(5.  42)  -4-2,o  ^  (1  —  w2)/4„  n  —  0, 1,  2,  •  ■  • 

and  the  An(Z2)  are  analytic  functions  of  Z2  which  are  dominated  by 
A/(R  —  Z2)n+ 2,  (i.  e.,  for  which  An(Z2)  «  A/(R — Z2)n+ 2  n  =  ' — 1,0,1, 
2.  •  •  •  holds)  A  and  R  being  suitably  chosen  positive  constants.  Then  the 
expression 

(5.  43)  if  =  Zx*  2  Bn^Z^Zi11,  s  =  1,  2, 

n= 0 

represents  a  solution  of  (5.  35)  which  is  defined  in 

(5.44)  E[|  Z1  |  <  (B—\  Z2 1)/(1  +  e), 

\Z,\  <1— ((— 1)*(1  —  44_2,  „)&  +  A)/(  1  - A-2, „)»] 

where 

(5.45)  e  =  max  [|(1 —  ( — l)s(l —  44-2,0)i  +  ^)/(l —  44-2l0)s  | , 

|  4/(1+  (— 1)*(1-4A*0)»)|] 

r.-i+  (—  l)s(i  —  4-2.o)5,  0<arg(i  —  A.Sj0)*O/2,  s-1,2. 

Here  B0(Z2)  is  an  arbitrary  function  for  which 

(5.46)  S0(Z2)  <<£(1  — tf-%)-2,  5>0, 

cmd  Bn(Z2)  are  to  be  determined  by  the  relations 

(5.47)  [1+  (—  1)*(1  —  4^_2,o)^]B1(Z2)  =  —  A^(Z2)Bo(Z2) 
n[n  -\-  ( — l)s(l  —  44._2j0)2]5n(^2) 

—  — \B"n-2{Z2)  -|-  2  4.v_2(Z2).Bn_v(Z2)],  n  =  2;  3,  •  • 

y=i 

We  should  like  to  add  here  a  remark  regarding  the  general  question  of 
analytic  continuation  of  a  solution  \p(Z,Z)  of  a  linear  partial  differential 
equation.  If  ^  is  given  in  two  different  domains,  say  Bx  and  B2  by  different 
representations,  say,  in  Bx  by  the  integral  operator  of  the  first  kind  in  the 
form 

(5.48) 

[see  (3. 11)]  and  in  B2  by  another  operator 

(5.49)  —  CE(Z,Z,t)f(iZ(l-t°))(dt/(l-t*)t) 

J  ^-1 
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[see  (3.4a)]  (not  necessarily  of  the  first  kind),  and  the  origin  lies  in 
Bi  Pi  B2}  then  the  problem  of  the  analytic  continuation  of  if/2  into  the  domain 

Bx  is  equivalent  to  the  determination  of  g  and  h  from  a  given  /=  2  araZ” 

__  n= o 

Setting  Z  =  0,  and  Z  =  0,  respectively,  in  the  relation 

(5.50)  p[?(2)]  +p[h(Z)\  =  riB(Z,2,i)/(iZ(l_i2))«/(i_^)5 

r=-i 

we  obtain  the  identities 


(5.51) 


where 
(5.  52) 


2  Zn 2  Tnv^&v  =  g (Z)  -{-  .R(0, Z) h (0), 

00 

aoo2^rn„<2>  —  R(0,Z)g{0)  +  h(Z) 


n~o  v~o 


Tn 


v(1)  =  f 1  En-v(1)  (1  —  t2)  v~*%~vdt, 

J  t=- 1 

j  t~~  i 


Tnn 

00 


E(Z,0,t)  =~2En(»(t)Z’>,  E(0,Z,t)  ='2Enm(t)Zn. 

71=0 


n=o 


[See  5,  page  310.]  It  is  thus  seen  that  the  analytic  continuation  of  if/2  into 
Bi  and  all  the  problems  arising  from  it,  such  as  the  determination  of  the 
singularities,  etc.,  are  reduced  to  similar  problems  in  the  theory  of  functions 
of  one  complex  variable  which  are  given  by  their  power  series  expansions. 


6.  The  determination  of  the  associate  function  in  terms  of  the  given 
streamf unction.  We  now  turn  to  the  problem  of  determining  the  associate 
function  f  in  terms  of  the  given  streamfunction.  In  the  case  of  the  integral 
operator  of  the  first  kind,  the  formulas  (3.  7),  (3.  10),  (3. 15),  (3.  16)  yield 
the  associate  g^iZ).  Therefore,  if  the  continuation  of  the  streamfunction 
(which  a  priori  is  given  in  the  real  plane)  to  complex  values  of  the  argu¬ 
ments  is  known,  then  these  formulas  immediately  yield  the  associate.  In  the 
case  of  general  integral  operators  (in  particular,  of  that  of  the  second  kind), 
if  the  streamfunction  is  given  in  the  form  of  a  power  series,  then  it  may  be 
shown  (see  [5],  p.  310)  that  the  associate  can  be  expressed  in  the  form  of  a 
power  series.  On  the  other  hand,  the  streamfunction  if/  is,  in  many  instances, 
given  in  some  different  form,  say  the  values  of  ip  and  dif//dM  are  given  on  a 
line  M  =  const.  If  these  quantities  are  analytic  functions  of  6 ,  then  from 
these  data  it  will  be  possible  to  determine  the  coefficients  of  the  series 
development  from  which  we  may  then  determine  the  associate  functions  in 
the  way  indicated  above.  However,  this  procedure  gives,,  only  the  function 
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element  of  f(Z),  so  that  the  function  is  in  general  determined  only  in  a 
sufficiently  small  neighborhood  of  the  interval  of  M  =  const,  along  which  the 
values  if/  =  Xi(0)  and  d\j//dM  =  x2(0)  are  given.  We  shall  show  in  the  present 
section  that,  in  the  case  of  certain  generating  functions  of  the  second  kind, 
a  formula  expressing  f  in  terms  of  the  values  </'  =  xiW  an^-  =  X2  W 

on  the  sonic  line,  M  —  1/  can  be  derived.  According  to  the  consideration 
of  5  an  integral  operator  of  the  second  kind  can  be  written  in  the  form  17 

(6.1)  if/(\,  6)  =  Im[  f  E(Z,Z,t)f(iZ(l-t2))(l-t*)-m 

J  c2 

=  (2  i)-1  r  [^/— ?)■*& 

J  c2 

where 

(6.2)  E  =  AXE^  +  HZ(1  —  t2)Y^A2E^2\  AK  complex  const. 

Here  EiK\  x  =  l,  2,  are  the  generating  functions  introduced  in  (5.5)  and 
(4.  1) ;  C2  is  a  simple  curve  in  the  complex  tf-plane  which  connects  t  = —  1 
with  t=l  and,  except  for  the  endpoints,  lies  outside  E[|  t  |  <  1].  Moreover, 
C2  has  to  be  chosen  in  such  a  manner  that  [JZ(1  —  t2)~\  lies  in  the  regularity 
domain  of  /  for  the  values  of  Z  under  consideration. 

We  assume  the  associate  function  to  be  of  the  form 

00 

(6.3)  /(£)  =  2  cvlv+1^,  cv  complex  const. 

I*=0 

which  is  suggested  by  previous  considerations.  Under  the  assumption  that 
the  (complex)  constants  Ax,  A2  satisfy  the  inequality 

(6.4)  Im[A2A1]  0, 

the  desired  inverse  formula  for  /  in  terms  of  xi  W  an(l  is  giyen  in 

Theorem  6.1.  Let  if/(\,0)  be  a  (real)  solution  of  the  compressi¬ 
bility  equation  (2.  8a)  which  is  defined  in  a  domain ,  say  E,  situated  in 
[31  |  x  |  <  |  6  | ,  X  <  0]  and  such  that  its  boundary  includes  an  interval 
[0O  <  e  <  0i]  of  the  transonic  line  X  =  0  not  containing  the  origin.  Let 

00 

(6.  5a)  lim  <p(\,  6)  =  Xi  W  —  2  avwP,  av(,)  real  const. 

X-*o-  *=0 

oo 

(6.5b)  lim  MKO)  “XsW  =3^(1  —  fe2)2/82  = 

X-»o-  •  *=° 

(see  (2.6))  t/ie  P2  associate  f  of  the  integral  operator  (6.1)  with 

generating  function  (6.  2)  and  At,  A2  satisfying  (6.  4)  is  given  by 

17  In  the  remainder  of  this  section  we  shall  omit  the  subscript  “  2  ”  in  Eit  and 
shall  write  simply  E. 
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(6.  6)  /({)  =  —  {  (—  2i01/Q/S1/^So2Im[A2A1]  }'[—  do  t~^Xl(a)dt 

J  c 

+  i  (— 1  )kdk  f  t-5/3Xk(<r)dt] 
fc= 1  C 

where 

(6.7)  .  cr  —  —  2t£  ( 1  —  t2) , 

the  constants  dv ,  v  =  0, 1,  2,  are 

(6.  8)  d0  =  —  (2/3)t»/*£oA2,  i  =  ~(25/3/3)^0/§1^1>  i2  „  —  i^S0A1} 
the  Sv  are  defined  in  (4.4). 

Remark  6. 1.  Developing  the  right-hand  side  of  (2.  6)  into  an  infinite 
series  and  inverting  it,  we  obtain  the  following  series  for  (6X/dM) 

(6.9)  dX/dM  =  [31/8 ( 1  —  &2) 2/3]  (—  X) i/3 

X  [l  +  31/3 (1  —  /t2) -4/3 ( (3/10)  -J-  J/t2  - —  (4/5)ft4)  ( —  A.)2/3  +  •  •  •] 

Therefore  it  is  legitimate  to  consider  lim  ( — A)1/3((ty/d A),  instead  of 

lira  m/dM).  x‘>0' 

A-»o- 

Proof.  A  formal  computation  yields  (see  (5.3),  (5.5),  (5.14),  (6.8)) 


(6. 10)  lim  E{Z,  Z,t)=—  4<-Vs^V« 

A— >o- 

and 

(6. 11)  lim  (—  \y/*Ex{Z,  Z,  t)  =  [dJ-V*  +  d0tsP(  1  —  <*)*/»]fl- >/“. 

A-»o- 

We  now  determine  the  limit  values  as  A-»  0“  for  the  right-hand  side  of  (6. 1) 
and  the  derivative  of  this  expression  multiplied  by  ( — A)1/3.  Since  by 

00 

assumption  these  limit  values  are  equal  to  expressions  2flv(1)0"  and 

00  V=0 

2  31/3(1  —  h2)2/*av{2)6v  respectively,  we  obtain  (by  equating  the  coefficients 

y=0 

of  the  0V)  the  following  system  of  equations 

(6. 12a)  av(1)  =  Im[( — 2i)~(m/6)<22/v(1)  cv], 

(6.12b)  uv(2>  =Im[(— 2t)-<y+1/e)(d1/v(1)  +d0Jv(2))cv],  v  —  1,  2,  ■  *  • 

where 


(6.13a)  Iv{1)  =  f  trVH\~t2y-xt*dt 

cs 

=  ±  J«-(4/3)»Ti(g-(4/8)Ti_1)  r(1/3)r(y  +  2/3)  n  =  0  ±1 
2  v  r(v  + 1) 


(6. 13b)  Ir(!>  —  f  t-5/3(l  — 


=  ±  _ 1) 


r(—  i/3)r(v  +  4/3) 
r(v  + 1) 
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We  note  that  the  last  term  on  the  right  hand  side  of  (6. 13b)  is  obtained 
as  follows: 

/v(2)  (3/2)  f  (1  - 72)^/3<7(7-2/3) 

c2 

=  —  (3/2)[*-2/3(l  —  3(v+l/3)  f  p!\l—t2)v-2!zdt. 

*•=-1  J  c2 

The  first  term  in  the  last  expression  vanishes  and  in  the  second  term,  as  in 
7v(1)?  we  replace  the  integration  curve  C2  by  the  segments  ( —  1,  0")  and 
(0%  1)  and  a  half-circle  around  the  origin. 

Eemark.  Since  the  curve  C2  need  only  satisfy  the  inequality 
|  t  |  >  |  2A/(A  +  i$)  |  *  (see  4. 11),  it  is  valid  to  replace  it  by  a  curve  con¬ 
sisting  of  the  segments  ( — 1,  0")  and  (0+,  1)  and  a  half -circle  around  the 
origin  provided  the  radius  of  the  latter  is  greater  than  ]  2 A/ (A  +  ^0  | 1  which 
approaches  zero  as  A-»0.  The  right  hand  sides  of  (6.13a)  and  (6.13b) 
are  obtained  by  making  the  substitution  t2  =  r  and  considering  the  integrals 
in  the  three-sheeted  r-plane.  In  the  following,  we  choose  the  sheet  for 
which  n  —  1. 

We  introduce  the  integrals 
(6.14a)  Jvw 


(6.14b)  Jv{2) 


and  note  that 

(6. 15) 

The  determinant  of  the  system  (6. 12)  does  not  vanish  for  any  v  since 
it  equals 

(6.16)  —  {i/2)lm[d0d2iw7&] 

which  would  vanish  only  if  Im[cZ0d2]  =  0.  Since,  by  (6.8),  Im[d0<f2] 
=  (2/3)$02  Im[A2Ai],  the  determinant  does  not  vanish,  because  of  condition 
(6.4). 


=  J  —  t2) 


rh  Je-(2/3)™^(e-(2/3),H  —  q) 


j  i—v) 


r(— i/3)r(v  + 1) 
r(v  +  s/3) 


=  ±  i e-(i/»)  »*•<(«- (*/*)*■<  —  i) 


r(l/3)r(y  +  i) 
r(v  +  4/3) 


lvwjvw  =  lyWjJ*)  =  —  3 
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fa  „  —  (—  2i)v+1/6  „  .  2 

1  ’  ImUA]  [~ *’  +  2  (- 

from  which  (6.6)  follows. 
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APPENDIX  I. 


The  derivation  of  formulae  (2.4).  In  this  section 
we  shall  at  first  determine  the  series  development 
of  -£(n)  -  — Z~Jl-  introduced  in  (2.2)  in  the  neigh- 

yz 

borhood  of  the  point  H  =  0.  According  to  (42)  of 
[6]  E  is  defined  by 


(iul)  fr  =  *  h(^(1))  =  o 

1 

( ll  2 

•q'1  } ~  [2/(k+l)]  being  the  speed  which  corresponds 
to  the  Mach  number  1. 

According  to  p,  7  of  (.6]  it  holds  for 


(A.  2)  T2  =  l-M2= 

5" 


I-  l/2(ict-l,.q2  > 

1-  l/2(k-l)q 

?  (k-1) 

(l-l/2(k-l>i2) 


Eliminating  v  between  both  equations  (A. 2)  we  get 

_1 _  _  _1 _ 

(A. 3)  ?  =  (~-)k_1  (1-  (— — -)?2)  ~  k-  1  . 

1 


IA.3)  P=  (—-V-1  (1-  (— — 

J  k+1  k  +  1 


=  (— )k_1  (1+ 

k+1 

Now,  from  (A.l)  we  have 

l  2 


k  T 


+  •  •  *) 


k+1  2 (k+1)" 

d.H  B  dH  d';q  =  dq,  _ 
dT 


d:i 


dT 


* 


dT 


=  jL  — L_  ,  prom  the  first  formula  of  (A. 2) 

2  dT 


.  ,  .  1  2 
we  obtain  g  : 


1  -  T _ 

(k+l)-  (k-l)  T2 


893, 

*  Therefore 


(A.4) 


dT 


-  2.  T.  +  2(k~l) . < 

1-  T2  (k+l)-(k-l)T! 


2_g T _ =  _  2  r _ T_  7  u 

(1-T^ (k+l)-(  k-l)  T2  k+  1  j  (l-T^(l-h2 


2 

(— ) 

k+l 


if  m2 

1+  _1_ 
k+l 


♦  + 


2(k+l)' 


1 

•  *  •  4 


-(— — > k- 1  rT*(Sij.)Ts.  t5 

k+l  L  k  +  1 


(k+l)' 


Therefore  (since  H=0  for  T=0) 

k  r_2 


(A. 5)  H=  -(  2  )  krr  E  +  (2k+i)rt  3JLTi 


k+l  L2  4(k+l)  6 (k+l) 

How  from  (A. 3)  follows j 


.5k+^T6  1 

k+T)2  J 


(A.  6) 


=  (JL-)k‘  1 

t2  k+  1; 


1,  2  1.5k  +2k+l  2 

-?  + - + - T  + 

T  k+l  ' 


(k+l) 


For  simplicity,  let 


(A.  7)  A= 


-  p'ihl 

o  L  -* 


k 

lirrT 


r,  $  ,  2  ,  k-  1 

B=  ~2“  ( - -) 

T  k+l 


_JL_ 

^(H)  'k  +  1 ' 


k-l 


Then 

/  2  /2k+  1  \4  3k2  +  2.5k  +  1  6 

(A.  8)  A=  T  +  ( - - -)T  + - r - — - T  + 


2k  +  2 


3(k+  l)2 


and 
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(A.9)  B  =  + 


1.2.  1.5k  +2k+l  „2 


k+1 


(k+1)' 


T  +  ... 


Nov/  (A. 8)  can  be  inverted  to  yield  as  a 

function  of  A  and  (A.9)  can  be  "turned  upside 

down"  to  give  i  as  a  Taylor  series  in  T 2.- 

There  are  thus  obtained  the  following  series* 

2k-l 
Fk+2 


(A. 10)  T2  -A  -  Ife^A2 


6k +6 


(A.U)  i  -  T2  ♦  *  I5  .  T6.  ...  . 

(k.l)4 


2  . 


Replacing  T  in  (A. 11)  by  (A. 10)  we  obtain 
(A. 12)  1(H)  -  i  =  A  r5k*5_)A2  + 

JV(k+fJ]  2/(k"1)  B  ^  Zk+Z 

+  k4+(43/6)k°+16k2+15.5k+(31.6)  3 

— 11  **’*"“ A  11  .  A  **■  •  «  • 

(k+1)4 

from  which  (2,4)  follows. 

The  derivation  of  formulae  (2,16)  and  (2,17)* 
fte  denote  in  the  following  by  u  the 
quantity 

(A. 13)  u  =  qW-  qs^)7  -  . 


It  follows  from  (A. 2)  that  in  the  neighborhood 
2  2 

of  u  «  0,  T  » 1  -  M  has  the  series  development 


Thr.s 

(A..  13)  A -'-3 


7 

>1)(W)4 


+ 
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+2  ^  3(k+l)^  lu  +  . .. 


Raising  each  side  to  2/3  power,  we  get 
12  5  2 


(A. 19) 


3j  (k+l)  6(-ft)3  = 


r  ~  _! 

u<L-2  5 

c 


(k+l)2(lc-5/2)u  +  ...j 


Inverting  (A. 19)  we  obtain 

JL  2  5 

(A. 20)  u  =  'q(l)-q  =  2  633(lc-l)  6(- ft) 


A'  1  i  ^  O 

+2  25"1(k+l)2(k-5«2)|^  63&(k+l)  6(-ft)3  2+... 
According  to  (47)  of  [6], 
N=-i(k+l)q4ji-i(ktl)q2  "  ’^^(k-Dq^T 

Expanding  this  expression  in  a  power  series  of 

u  =  2**^- A  (k+l  we  get 

J9  5  ^  3p  p  ,1 

(A. 21)  H=-2  *(k+l)"\  ^il-|i(k-l ) (k+l)2  + 

1  ^  L  In  n 

2^ (13/8) (k+l)2  u  +  ... J 

Raising  teth  sides  to  the  (~2/s) ‘ power,  we 


obtain 


£  3  ^5  .*  p 

3-  o^Vwn wSj.,,.  A 


(A. 22)  N  =  2‘,(k+l)  |u+j~(k-l)(k+l)' 

H  22(k+l)2!u&  +  "4 


897. 


Replacing  u  by  the  expansion  (A. 20),  we  have 

-1  JL  J:  f£  x  J5  £t  J. 

2.  2(k+l)  bH  3  -j332~^(k+l)  3(-A)3J+5"12  2(k--2.S) 

X  332.  6(k+l)  S(-r^)3  *  ~(k-l)(k+l)2  + 

o 


g  22(k+l)  j  b32  6(k+l)  6(-A)3J  +  ...J 


_£  5  5  ?  _1  _5  Zf  ; 

(A. 23)  II  3=  22"(k+l)5  332  6(k+l)  6(-/l)3>b.+(k+l) 


rc^)«-iair44 


5  2 

6/  a\3 


|2  °S*(k+l)  °(^r\  +  .e. 


Raising  both  sides  to  the  (-3/2)  power,  we  obtain 


1  ~  ,  3 ,  .  2  ~T 

-I  12---2  ?. 


{A.  24)  *4, 


2.  633(k+l)  6(-A)3  • 


According  to  (A. 24)  and/ 70 )  of  [6],  we  obtain 


(A. 25)  P--(H24^=  —  0(;»  3  + 


A_?(-^  -  Ao  +  A2(-A)3  +  ... 


Remark.  In  particular,  for  k  =  1*401  we  have 


898 


(A. 25a) 


t?  =  V'"  A  u 

fej  A=»5  • 


where  A  =  ,45,  A  =  0,  A  -  .41C 73,  A  =  .32792, 
*0  -4  U 

A0  =  11-676.  A.*  .05392,  A -..0.1759,  A  -.09523, 

2  4  o  o 

Alr~  -,00149,  etc. 


3.  The  derivation  of  formula  (4,3),  (4,4) 

Expanding  the  right  hand  side  of  (?,6j  in  a  powe: 
series,  we  obtain 

(A.2,<S)->,^(l'-h2)T34(l-hi)T5-'-  -(l-dtr7  +  ...  , 


(A, 27)  f.,3^/(l..hZ)]-^ 


Then?  ^ 

(A. 20)  =T+  iCl+h2)!**  ..|g.  {28+llh2+18b4.)T6  + 

Inverting  thij  series  we  ojtala.. 

(A. 29)  T  =  ^3-  I(l+h2)^b12)A  ~~  (3+21h^5h4/>7) 

±(  I 

(A, 30)  t'“  ^  3-jl+  ici+h2)^3)2  + 


+  Ag(4*-17h2+4h  4)(yx5)4  +  ...  j- 


or  finally 


1  _ir 

2  el.  i  /, 2 


(A. 31)  T  «yw_  un+  Y^(l+h  )(^13)  + 

^  I  ? 

T40V(9-82h  +  »b)CMn3)4+  —  < 
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Now 


(A. 32)  £[+  i(k-l)M^ 


1 

2(k-l) 


r  i 


2/(lc+l  )i  ^FT) 


l-h2T2J~ 


2/(k+l) 


2(k-mi+  l(k-1  )“1h2T2*  4(k-l)~2(2k-l)h4T4 


0  2 

1  ,  ^  -3 _  _ 6 


*  ■JgCk-l)"  (2k-l)(4k-3)h 


Replacing  T  in  (A. 32)  by  its  power  series  expansion 
(A. 29)  we  obtain 


(A. 33)  jl  +  i(k-l)M2J ~2(k-l)=|”g/(k+i 


1 

2TCT) 


ft.  1 

At  i(k+l  r^uf-  ^(k+1  )_2-(6k*5 ) 


r\ 

-giM(k+l)"S'(584k2  +  122k  -  579)  L  + 

/  sJ 


If  we  multiply  series  (A,31)  and  (A. 33)  together, 
we  get  the  following  expansion  for  11(2^,)  as  a  function 
°f  =  -3^/(l-h2)  =  -  §(k+l)ft  , 

(A..34)  H(2  /\)  =  2(k+l)2^*Iyi  6/l  + 


W<W)'V*5)ft!  -IJoo0"1)'2 


(64k‘ 


'+  70k  +  75)^13  +  ... j 


900. 


Derivation  of  formula  (6.9  ).  Differentiating 
(A, 26)  with  respect  to  T,  we  have 

(A. 35)  -(d//dT)=(l-h2)T2*(l-h4)T4-*!(l-h6)T6  + 

2  2 

Since  T  =  1-M  and  therefore 


(A.36)-(dT/dia)®T_1M=T~1  ^1--|t2- il4  +  ...j| 

we  obtain  by  multiplying  (A. 35)  by  (A. 36) 

(A. 37)  (dydM)"(l-h2)T+Og-  +  -|h2  -  h4)TS 

^  ,3  1.2  1,6  .  6,m5 

+  (8  +  8h  +  2h  "h  >T  +  “• 


Replacing  in  (A. 37)  T  by  its  series  development 
(A. 29)  we  find 


li 


(A. 38)  (dVdU)-(l-h ?)j\°  +  +  Ih2“ 

+  + 

5.  Formula  for  and  *  A  straightforward 

^  A  *0  e 

derivation  of  (4.15)  with  respect  to  A  and  0,  res- 
pectively,  yields 

(A  S9)^l-  P 

7>'h  1 


l 
3x2 


(-2ft)0  (-t  (7\+  10)) 

(p ,(»  ^  i 

..  3  6  %*  i® 


"T 

-t2(?v 


_ 'X* 

( ?v+  is/ 


_1 

*\\  3. 


+  2i(-27Q°Q  A  2  1 


t2(A+  i9)2 


F  f  VnS’fzS’TT* 


6’3-’3’  ,M 

^  .t  (A+  10), 


v+  19)')' 
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(t  (A*ie)) 


'  5  W) 

p Aa)-  *  — ? . . 

6  (-t2(Vie) 


p/S  4  5  -2  A  Y  2i(^2A)9 

_t2(^+ie)y  t^+ie)2 


jfs  4  5  -2  A 

6  3  3  -t  (  A+iQ) 


where 


P3(A)=-jj  Sj+j  S2(-2^3+  ..^f+sx(-2^)+  ..  J 

P4(/0=-p  +|si(-2^)5+  ...1  ^/|l+31(-2A)3+...J 


(A.40^~ 


= _ mL  4  .  1  LpA  2  1  -2  A 

c*  (-t2)*(>ie)3  (,2[  6  '6"3J3’-t2(?,+ie). 


_L4__.F7i  2  i  ±% _ A  - 

t2(/^:.G)  '6’3’3’  -t2(>\-.'-ie)/ 


(-2?) 


-B2  -  6  -I 


•t2)°  (/N+iG) 


— — -V 

+iflA  -6  3  3  -^(A+iS)^ 


- F’/5  4  5  -2  A  \  I  ( 

( „t2 )#(A+T^7  6 ’ 3 ' 3 '  -t2 (>»+  ie)'J  ) 


ERRATA 


u  n 

pg,  860,  1.  4  from  bot.t  i  should  be  deleted, 
pg*  862,  formula  (3.2)  should  read 

tt 

"  ^[(l-t2)E(Z,Z,t)/zt]/^  Z 

n  it 

PK.  872,  formula  (4.9*);  the  exponent  -2/3 

Tt  tt 

should  be  changed  to  2/3  . 

tt  tt 

pg.  872,  formula  (4,11)  \  (  should  be  added 

n  11 

after  dt/  . 

pg.  882,  formula  (5.29):  the  symbol  before  S 

it  if  '  y 

is  =  ;  an  s  should  be  deleted, 

pg.  892,  line  5  from  bot.  should  read 

tt  11 

Eliminating  q  between  both  .  .  , 

pg.  892,  line  2  from  bot,  should  read 
Now,  from  (A.l)  we  have  dT  q  ^ 
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